This is the analysis script to prepare the data to run analyses for the project “Anxiety, depression and trauma symptom change during the COVID-19 pandemic: retrospective versus objective assessment” - Young et al (2020)

Script written by K L Purves, K N Thompson, C Huebel and M Davies. Email: kirstin.purves@kcl.ac.uk, katie.thompson@kcl.ac.uk, christopher.1.huebel@kcl.ac.uk, molly.davies@kcl.ac.uk

#Set up

Clear global environment

Retrieve the current date to use it for file endings to not overwrite files when one exports files

#Packages Install packages (if they are not available in your version of R)

Load packages

R functions

Half-violin plot

Colour palettes

Define colours for plotting this are the standard coping colours

Choose in this chunk which palette to use 04.07.2020 - Default to 2 colour COPING palette

ggplot theme

Set up ggplot theme for the plots

Function for split violin plot

Data import

Source data file paths

Read in data

[1] 31884   641

  [1] "Sample"                                                           "ID"                                                              
  [3] "startDate.prepandemic"                                            "Age_uncleaned"                                                   
  [5] "Gender_unc"                                                       "Sex"                                                             
  [7] "EduYrs"                                                           "demographics.college_or_university_degree"                       
  [9] "demographics.a_levelsas_levels_or_equivalent"                     "demographics.o_levelsgcses_or_equivalent"                        
 [11] "demographics.cses_or_equivalent"                                  "demographics.nvq_or_hnd_or_hnc_or_equivalent"                    
 [13] "Ethnicity_unc"                                                    "Gender_numeric"                                                  
 [15] "Sex_numeric"                                                      "demographics.college_or_university_degree_numeric"               
 [17] "demographics.a_levelsas_levels_or_equivalent_numeric"             "demographics.o_levelsgcses_or_equivalent_numeric"                
 [19] "demographics.cses_or_equivalent_numeric"                          "demographics.nvq_or_hnd_or_hnc_or_equivalent_numeric"            
 [21] "Ethnicity_numeric"                                                "Gender"                                                          
 [23] "Ethnicity"                                                        "Age"                                                             
 [25] "age_category_numeric"                                             "age_category"                                                    
 [27] "highest_education_prepan_numeric"                                 "highest_education_prepan"                                        
 [29] "employment.phd"                                                   "employment.masters_degree_or_equivalent"                         
 [31] "employment.postgraduate_degree_or_equivalent"                     "employment.bachelors_degree_or_equivalent"                       
 [33] "employment.a_levelsas_levels_or_equivalent"                       "employment.o_levelsgcses_or_equivalent"                          
 [35] "employment.cses_or_equivalent"                                    "employment.nvq_or_hnd_or_hnc_or_equivalent"                      
 [37] "employment.government_work_key_workers"                           "employment.unemployed"                                           
 [39] "employment.zerohours_contract"                                    "employment.stayathome_parent_or_carer"                           
 [41] "employment.retired"                                               "employment.student_gcse_or_a_level"                              
 [43] "employment.student_university"                                    "employment.became_unemployed"                                    
 [45] "employment.furloughed_or_paid_leave_government_funded"            "employment.furloughed_or_paid_leave_company_funded"              
 [47] "employment.paid_leave_furloughed"                                 "employment.fulltime_employed"                                    
 [49] "employment.parttime_employed"                                     "employment.selfemployed"                                         
 [51] "employment.contract_or_freelance_work"                            "employment.small_business_owner"                                 
 [53] "employment.receiving_state_income"                                "employment.my_employment_status_has_not_changed"                 
 [55] "employment.reduction_in_hours"                                    "employment.reduction_in_salary"                                  
 [57] "employment.benefits_increased"                                    "employment.benefits_decreased"                                   
 [59] "employment.became_employed"                                       "employment.increased_hours"                                      
 [61] "employment.increased_salary"                                      "employment.phd_numeric"                                          
 [63] "employment.masters_degree_or_equivalent_numeric"                  "employment.postgraduate_degree_or_equivalent_numeric"            
 [65] "employment.bachelors_degree_or_equivalent_numeric"                "employment.a_levelsas_levels_or_equivalent_numeric"              
 [67] "employment.o_levelsgcses_or_equivalent_numeric"                   "employment.cses_or_equivalent_numeric"                           
 [69] "employment.nvq_or_hnd_or_hnc_or_equivalent_numeric"               "employment.government_work_key_workers_numeric"                  
 [71] "employment.unemployed_numeric"                                    "employment.zerohours_contract_numeric"                           
 [73] "employment.stayathome_parent_or_carer_numeric"                    "employment.retired_numeric"                                      
 [75] "employment.student_gcse_or_a_level_numeric"                       "employment.student_university_numeric"                           
 [77] "employment.became_unemployed_numeric"                             "employment.furloughed_or_paid_leave_government_funded_numeric"   
 [79] "employment.furloughed_or_paid_leave_company_funded_numeric"       "employment.paid_leave_furloughed_numeric"                        
 [81] "employment.fulltime_employed_numeric"                             "employment.parttime_employed_numeric"                            
 [83] "employment.selfemployed_numeric"                                  "employment.contract_or_freelance_work_numeric"                   
 [85] "employment.small_business_owner_numeric"                          "employment.receiving_state_income_numeric"                       
 [87] "employment.my_employment_status_has_not_changed_numeric"          "employment.reduction_in_hours_numeric"                           
 [89] "employment.reduction_in_salary_numeric"                           "employment.benefits_increased_numeric"                           
 [91] "employment.benefits_decreased_numeric"                            "employment.became_employed_numeric"                              
 [93] "employment.increased_hours_numeric"                               "employment.increased_salary_numeric"                             
 [95] "highest_education_numeric"                                        "highest_education"                                               
 [97] "employment_prior_covid_numeric"                                   "employment_prior_covid"                                          
 [99] "employment_change_numeric"                                        "employment_change"                                               
[101] "key_worker_numeric"                                               "key_worker"                                                      
[103] "anxiety_disorders_numeric"                                        "anxiety_disorders"                                               
[105] "mhd.specific_phobia_e.g._phobia_of_flying_numeric"                "mhd.specific_phobia_e.g._phobia_of_flying"                       
[107] "mhd.social_anxiety_or_social_phobia_numeric"                      "mhd.social_anxiety_or_social_phobia"                             
[109] "agoraphobia_panic_disorder_numeric"                               "agoraphobia_panic_disorder"                                      
[111] "depressive_disorders_numeric"                                     "depressive_disorders"                                            
[113] "eating_disorders_numeric"                                         "eating_disorders"                                                
[115] "obsessive_compulsive_disorders_numeric"                           "obsessive_compulsive_disorders"                                  
[117] "psychotic_disorders_numeric"                                      "psychotic_disorders"                                             
[119] "mhd.mania_hypomania_bipolar_or_manicdepression_numeric"           "mhd.mania_hypomania_bipolar_or_manicdepression"                  
[121] "mhd.posttraumatic_stress_disorder_ptsd_numeric"                   "mhd.posttraumatic_stress_disorder_ptsd"                          
[123] "autism_spectrum_disorder_numeric"                                 "autism_spectrum_disorder"                                        
[125] "mhd.attention_deficit_hyperactivity_disorder_numeric"             "mhd.attention_deficit_hyperactivity_disorder"                    
[127] "mhd.personality_disorder_numeric"                                 "mhd.personality_disorder"                                        
[129] "personality_cluster_a_numeric"                                    "personality_cluster_a"                                           
[131] "personality_cluster_b_numeric"                                    "personality_cluster_b"                                           
[133] "personality_cluster_c_numeric"                                    "personality_cluster_c"                                           
[135] "depression_and_anxiety_numeric"                                   "depression_and_anxiety"                                          
[137] "restricting_vs_binge_eating_eating_disorders_numeric"             "restricting_vs_binge_eating_eating_disorders"                    
[139] "bipolar_and_schizophrenia_numeric"                                "bipolar_and_schizophrenia"                                       
[141] "control_numeric"                                                  "control"                                                         
[143] "disorders_total_count"                                            "comorbidity_total_count_numeric"                                 
[145] "comorbidity_total_count_factor_numeric"                           "comorbidity_total_count_factor"                                  
[147] "mhd.none_of_the_above"                                            "mhd.none_of_the_above.1"                                         
[149] "gad.feeling_nervous_anxious_or_on_edge_prepan"                    "gad.control_worrying_stop_prepan"                                
[151] "gad.worrying_too_much_about_different_things_prepan"              "gad.trouble_relaxing_prepan"                                     
[153] "gad.sit_restless_hard_prepan"                                     "gad.becoming_easily_annoyed_or_irritable_prepan"                 
[155] "gad.awful_feeling_afraid_happen_prepan"                           "gad.feeling_nervous_anxious_or_on_edge_prepan_numeric"           
[157] "gad.control_worrying_stop_prepan_numeric"                         "gad.worrying_too_much_about_different_things_prepan_numeric"     
[159] "gad.trouble_relaxing_prepan_numeric"                              "gad.sit_restless_hard_prepan_numeric"                            
[161] "gad.becoming_easily_annoyed_or_irritable_prepan_numeric"          "gad.awful_feeling_afraid_happen_prepan_numeric"                  
[163] "gad.sum_score_prepan_raw"                                         "na_per_person_gad_prepan"                                        
[165] "gad.sum_score_prepan"                                             "pcl.stressful_experience_repeated_images_prepan"                 
[167] "pcl.stressful_experience_upset_reminded_prepan"                   "pcl.stressful_situation_avoiding_activities_prepan"              
[169] "pcl.cut_people_feeling_distant_prepan"                            "pcl.feeling_irritable_or_having_angry_outbursts_prepan"          
[171] "pcl.difficulty_concentrating_prepan"                              "pcl.stressful_experience_repeated_images_prepan_numeric"         
[173] "pcl.stressful_experience_upset_reminded_prepan_numeric"           "pcl.stressful_situation_avoiding_activities_prepan_numeric"      
[175] "pcl.cut_people_feeling_distant_prepan_numeric"                    "pcl.feeling_irritable_or_having_angry_outbursts_prepan_numeric"  
[177] "pcl.difficulty_concentrating_prepan_numeric"                      "pcl.sum_score_prepan_raw"                                        
[179] "na_per_person_pcl_prepan"                                         "pcl.sum_score_prepan"                                            
[181] "phq.little_interest_or_pleasure_in_doing_things_prepan"           "phq.feeling_down_depressed_or_hopeless_prepan"                   
[183] "phq.staying_asleep_sleeping_trouble_prepan"                       "phq.feeling_tired_or_having_little_energy_prepan"                
[185] "phq.poor_appetite_or_overeating_prepan"                           "phq.feeling_bad_failure_family_prepan"                           
[187] "phq.trouble_concentrating_reading_newspaper_prepan"               "phq.moving_fidgety_noticed_opposite_prepan"                      
[189] "phq.dead_hurting_thoughts_prepan"                                 "phq.little_interest_or_pleasure_in_doing_things_prepan_numeric"  
[191] "phq.feeling_down_depressed_or_hopeless_prepan_numeric"            "phq.staying_asleep_sleeping_trouble_prepan_numeric"              
[193] "phq.feeling_tired_or_having_little_energy_prepan_numeric"         "phq.poor_appetite_or_overeating_prepan_numeric"                  
[195] "phq.feeling_bad_failure_family_prepan_numeric"                    "phq.trouble_concentrating_reading_newspaper_prepan_numeric"      
[197] "phq.moving_fidgety_noticed_opposite_prepan_numeric"               "phq.dead_hurting_thoughts_prepan_numeric"                        
[199] "phq.sum_score_prepan_raw"                                         "na_per_person_phq_prepan"                                        
[201] "phq.missing_only_suicide_item_prepan"                             "phq.sum_score_prepan"                                            
[203] "phq.sum_score_8items_prepan"                                      "ocir.i_repeatedly_check_doors_windows_drawers_etc_prepan"        
[205] "ocir.arranged_things_change_upset_prepan"                         "ocir.repeat_feel_numbers_prepan"                                 
[207] "ocir.clean_simply_wash_feel_prepan"                               "ocir.mind_upset_unpleasant_thoughts_prepan"                      
[209] "ocir.afraid_avoid_throwing_things_prepan"                         "ocir.light_switches_water_taps_prepan"                           
[211] "ocir.arranged_things_prepan"                                      "ocir.good_feel_bad_numbers_prepan"                               
[213] "ocir.hands_longer_wash_prepan"                                    "ocir.saved_things_prepan"                                        
[215] "ocir.frequently_difficulty_rid_nasty_prepan"                      "ocir.i_check_things_more_often_than_necessary_prepan"            
[217] "ocir.arranged_properly_objects_upset_prepan"                      "ocir.feel_compelled_count_things_prepan"                         
[219] "ocir.strangers_touch_touched_difficult_prepan"                    "ocir.control_difficult_find_thoughts_prepan"                     
[221] "ocir.i_collect_things_i_dont_need_prepan"                         "ocir.i_repeatedly_check_doors_windows_drawers_etc_prepan_numeric"
[223] "ocir.arranged_things_change_upset_prepan_numeric"                 "ocir.repeat_feel_numbers_prepan_numeric"                         
[225] "ocir.clean_simply_wash_feel_prepan_numeric"                       "ocir.mind_upset_unpleasant_thoughts_prepan_numeric"              
[227] "ocir.afraid_avoid_throwing_things_prepan_numeric"                 "ocir.light_switches_water_taps_prepan_numeric"                   
[229] "ocir.arranged_things_prepan_numeric"                              "ocir.good_feel_bad_numbers_prepan_numeric"                       
[231] "ocir.hands_longer_wash_prepan_numeric"                            "ocir.saved_things_prepan_numeric"                                
[233] "ocir.frequently_difficulty_rid_nasty_prepan_numeric"              "ocir.i_check_things_more_often_than_necessary_prepan_numeric"    
[235] "ocir.arranged_properly_objects_upset_prepan_numeric"              "ocir.feel_compelled_count_things_prepan_numeric"                 
[237] "ocir.strangers_touch_touched_difficult_prepan_numeric"            "ocir.control_difficult_find_thoughts_prepan_numeric"             
[239] "ocir.i_collect_things_i_dont_need_prepan_numeric"                 "ocir.sum_score_prepan_raw"                                       
[241] "na_per_person_ocir_prepan"                                        "ocir.sum_score_prepan"                                           
[243] "startDate.coping"                                                 "gad.problems_made_difficult_care"                                
[245] "gad.pandemic_felt_feelings"                                       "gad.feeling_nervous_anxious_or_on_edge_retro"                    
[247] "gad.control_worrying_stop_retro"                                  "gad.worrying_too_much_about_different_things_retro"              
[249] "gad.trouble_relaxing_retro"                                       "gad.sit_restless_hard_retro"                                     
[251] "gad.becoming_easily_annoyed_or_irritable_retro"                   "gad.awful_feeling_afraid_happen_retro"                           
[253] "gad.feeling_nervous_anxious_or_on_edge_base"                      "gad.control_worrying_stop_base"                                  
[255] "gad.worrying_too_much_about_different_things_base"                "gad.trouble_relaxing_base"                                       
[257] "gad.sit_restless_hard_base"                                       "gad.becoming_easily_annoyed_or_irritable_base"                   
[259] "gad.awful_feeling_afraid_happen_base"                             "gad.problems_made_difficult_care_numeric"                        
[261] "gad.pandemic_felt_feelings_numeric"                               "gad.feeling_nervous_anxious_or_on_edge_retro_numeric"            
[263] "gad.control_worrying_stop_retro_numeric"                          "gad.worrying_too_much_about_different_things_retro_numeric"      
[265] "gad.trouble_relaxing_retro_numeric"                               "gad.sit_restless_hard_retro_numeric"                             
[267] "gad.becoming_easily_annoyed_or_irritable_retro_numeric"           "gad.awful_feeling_afraid_happen_retro_numeric"                   
[269] "gad.feeling_nervous_anxious_or_on_edge_base_numeric"              "gad.control_worrying_stop_base_numeric"                          
[271] "gad.worrying_too_much_about_different_things_base_numeric"        "gad.trouble_relaxing_base_numeric"                               
[273] "gad.sit_restless_hard_base_numeric"                               "gad.becoming_easily_annoyed_or_irritable_base_numeric"           
[275] "gad.awful_feeling_afraid_happen_base_numeric"                     "gad.sum_score_retro_raw"                                         
[277] "gad.sum_score_base_raw"                                           "na_per_person_gad_base"                                          
[279] "na_per_person_gad_retro"                                          "gad.sum_score_base"                                              
[281] "gad.sum_score_retro"                                              "pcl.stressful_experience_repeated_images_base"                   
[283] "pcl.stressful_experience_upset_reminded_base"                     "pcl.stressful_situation_avoiding_activities_base"                
[285] "pcl.cut_people_feeling_distant_base"                              "pcl.feeling_irritable_or_having_angry_outbursts_base"            
[287] "pcl.difficulty_concentrating_base"                                "pcl.stressful_experience_repeated_images_retro"                  
[289] "pcl.stressful_experience_upset_reminded_retro"                    "pcl.stressful_situation_avoiding_activities_retro"               
[291] "pcl.cut_people_feeling_distant_retro"                             "pcl.feeling_irritable_or_having_angry_outbursts_retro"           
[293] "pcl.difficulty_concentrating_retro"                               "pcl.stressful_experience_repeated_images_base_numeric"           
[295] "pcl.stressful_experience_upset_reminded_base_numeric"             "pcl.stressful_situation_avoiding_activities_base_numeric"        
[297] "pcl.cut_people_feeling_distant_base_numeric"                      "pcl.feeling_irritable_or_having_angry_outbursts_base_numeric"    
[299] "pcl.difficulty_concentrating_base_numeric"                        "pcl.stressful_experience_repeated_images_retro_numeric"          
[301] "pcl.stressful_experience_upset_reminded_retro_numeric"            "pcl.stressful_situation_avoiding_activities_retro_numeric"       
[303] "pcl.cut_people_feeling_distant_retro_numeric"                     "pcl.feeling_irritable_or_having_angry_outbursts_retro_numeric"   
[305] "pcl.difficulty_concentrating_retro_numeric"                       "pcl.sum_score_retro_raw"                                         
[307] "pcl.sum_score_base_raw"                                           "na_per_person_pcl_base"                                          
[309] "na_per_person_pcl_retro"                                          "pcl.sum_score_base"                                              
[311] "pcl.sum_score_retro"                                              "phq.problems_made_difficult_care"                                
[313] "phq.pandemic_felt_feelings"                                       "phq.little_interest_or_pleasure_in_doing_things_retro"           
[315] "phq.feeling_down_depressed_or_hopeless_retro"                     "phq.staying_asleep_sleeping_trouble_retro"                       
[317] "phq.feeling_tired_or_having_little_energy_retro"                  "phq.poor_appetite_or_overeating_retro"                           
[319] "phq.feeling_bad_failure_family_retro"                             "phq.trouble_concentrating_reading_newspaper_retro"               
[321] "phq.moving_fidgety_noticed_opposite_retro"                        "phq.dead_hurting_thoughts_retro"                                 
[323] "phq.little_interest_or_pleasure_in_doing_things_base"             "phq.feeling_down_depressed_or_hopeless_base"                     
[325] "phq.staying_asleep_sleeping_trouble_base"                         "phq.feeling_tired_or_having_little_energy_base"                  
[327] "phq.poor_appetite_or_overeating_base"                             "phq.feeling_bad_failure_family_base"                             
[329] "phq.trouble_concentrating_reading_newspaper_base"                 "phq.moving_fidgety_noticed_opposite_base"                        
[331] "phq.dead_hurting_thoughts_base"                                   "phq.problems_made_difficult_care_numeric"                        
[333] "phq.pandemic_felt_feelings_numeric"                               "phq.little_interest_or_pleasure_in_doing_things_retro_numeric"   
[335] "phq.feeling_down_depressed_or_hopeless_retro_numeric"             "phq.staying_asleep_sleeping_trouble_retro_numeric"               
[337] "phq.feeling_tired_or_having_little_energy_retro_numeric"          "phq.poor_appetite_or_overeating_retro_numeric"                   
[339] "phq.feeling_bad_failure_family_retro_numeric"                     "phq.trouble_concentrating_reading_newspaper_retro_numeric"       
[341] "phq.moving_fidgety_noticed_opposite_retro_numeric"                "phq.dead_hurting_thoughts_retro_numeric"                         
[343] "phq.little_interest_or_pleasure_in_doing_things_base_numeric"     "phq.feeling_down_depressed_or_hopeless_base_numeric"             
[345] "phq.staying_asleep_sleeping_trouble_base_numeric"                 "phq.feeling_tired_or_having_little_energy_base_numeric"          
[347] "phq.poor_appetite_or_overeating_base_numeric"                     "phq.feeling_bad_failure_family_base_numeric"                     
[349] "phq.trouble_concentrating_reading_newspaper_base_numeric"         "phq.moving_fidgety_noticed_opposite_base_numeric"                
[351] "phq.dead_hurting_thoughts_base_numeric"                           "phq.sum_score_retro_raw"                                         
[353] "phq.sum_score_base_raw"                                           "na_per_person_phq_retro"                                         
[355] "phq.missing_only_suicide_item_retro"                              "na_per_person_phq_base"                                          
[357] "phq.missing_only_suicide_item_base"                               "phq.sum_score_base"                                              
[359] "phq.sum_score_8items_base"                                        "phq.sum_score_retro"                                             
[361] "phq.sum_score_8items_retro"                                       "ocir.felt_pandemic_feelings"                                     
[363] "ocir.i_repeatedly_check_doors_windows_drawers_etc_retro"          "ocir.arranged_things_change_upset_retro"                         
[365] "ocir.repeat_feel_numbers_retro"                                   "ocir.clean_simply_wash_feel_retro"                               
[367] "ocir.mind_upset_unpleasant_thoughts_retro"                        "ocir.afraid_avoid_throwing_things_retro"                         
[369] "ocir.light_switches_water_taps_retro"                             "ocir.arranged_things_retro"                                      
[371] "ocir.good_feel_bad_numbers_retro"                                 "ocir.hands_longer_wash_retro"                                    
[373] "ocir.saved_things_retro"                                          "ocir.frequently_difficulty_rid_nasty_retro"                      
[375] "ocir.i_check_things_more_often_than_necessary_retro"              "ocir.arranged_properly_objects_upset_retro"                      
[377] "ocir.feel_compelled_count_things_retro"                           "ocir.strangers_touch_touched_difficult_retro"                    
[379] "ocir.control_difficult_find_thoughts_retro"                       "ocir.i_collect_things_i_dont_need_retro"                         
[381] "ocir.i_repeatedly_check_doors_windows_drawers_etc_base"           "ocir.arranged_things_change_upset_base"                          
[383] "ocir.repeat_feel_numbers_base"                                    "ocir.clean_simply_wash_feel_base"                                
[385] "ocir.mind_upset_unpleasant_thoughts_base"                         "ocir.afraid_avoid_throwing_things_base"                          
[387] "ocir.light_switches_water_taps_base"                              "ocir.arranged_things_base"                                       
[389] "ocir.good_feel_bad_numbers_base"                                  "ocir.hands_longer_wash_base"                                     
[391] "ocir.saved_things_base"                                           "ocir.frequently_difficulty_rid_nasty_base"                       
[393] "ocir.i_check_things_more_often_than_necessary_base"               "ocir.arranged_properly_objects_upset_base"                       
[395] "ocir.feel_compelled_count_things_base"                            "ocir.strangers_touch_touched_difficult_base"                     
[397] "ocir.control_difficult_find_thoughts_base"                        "ocir.i_collect_things_i_dont_need_base"                          
[399] "ocir.felt_pandemic_feelings_numeric"                              "ocir.i_repeatedly_check_doors_windows_drawers_etc_retro_numeric" 
[401] "ocir.arranged_things_change_upset_retro_numeric"                  "ocir.repeat_feel_numbers_retro_numeric"                          
[403] "ocir.clean_simply_wash_feel_retro_numeric"                        "ocir.mind_upset_unpleasant_thoughts_retro_numeric"               
[405] "ocir.afraid_avoid_throwing_things_retro_numeric"                  "ocir.light_switches_water_taps_retro_numeric"                    
[407] "ocir.arranged_things_retro_numeric"                               "ocir.good_feel_bad_numbers_retro_numeric"                        
[409] "ocir.hands_longer_wash_retro_numeric"                             "ocir.saved_things_retro_numeric"                                 
[411] "ocir.frequently_difficulty_rid_nasty_retro_numeric"               "ocir.i_check_things_more_often_than_necessary_retro_numeric"     
[413] "ocir.arranged_properly_objects_upset_retro_numeric"               "ocir.feel_compelled_count_things_retro_numeric"                  
[415] "ocir.strangers_touch_touched_difficult_retro_numeric"             "ocir.control_difficult_find_thoughts_retro_numeric"              
[417] "ocir.i_collect_things_i_dont_need_retro_numeric"                  "ocir.i_repeatedly_check_doors_windows_drawers_etc_base_numeric"  
[419] "ocir.arranged_things_change_upset_base_numeric"                   "ocir.repeat_feel_numbers_base_numeric"                           
[421] "ocir.clean_simply_wash_feel_base_numeric"                         "ocir.mind_upset_unpleasant_thoughts_base_numeric"                
[423] "ocir.afraid_avoid_throwing_things_base_numeric"                   "ocir.light_switches_water_taps_base_numeric"                     
[425] "ocir.arranged_things_base_numeric"                                "ocir.good_feel_bad_numbers_base_numeric"                         
[427] "ocir.hands_longer_wash_base_numeric"                              "ocir.saved_things_base_numeric"                                  
[429] "ocir.frequently_difficulty_rid_nasty_base_numeric"                "ocir.i_check_things_more_often_than_necessary_base_numeric"      
[431] "ocir.arranged_properly_objects_upset_base_numeric"                "ocir.feel_compelled_count_things_base_numeric"                   
[433] "ocir.strangers_touch_touched_difficult_base_numeric"              "ocir.control_difficult_find_thoughts_base_numeric"               
[435] "ocir.i_collect_things_i_dont_need_base_numeric"                   "ocir.sum_score_retro_raw"                                        
[437] "ocir.sum_score_base_raw"                                          "na_per_person_ocir_retro"                                        
[439] "na_per_person_ocir_base"                                          "ocir.sum_score_base"                                             
[441] "ocir.sum_score_retro"                                             "alsfrs.limb_weakness_noticed_past"                               
[443] "alsfrs.normal"                                                    "alsfrs.slow_or_sloppy_all_words_are_legible"                     
[445] "alsfrs.not_all_words_are_legible"                                 "alsfrs.able_to_grip_pen"                                         
[447] "alsfrs.unable_to_grip_pen"                                        "alsfrs.have_you_noticed_any_difficulty_walking"                  
[449] "alsfrs.do_you_have_any_difficulty_climbing_stairs"                "alsfrs.limb_weakness_noticed_past_numeric"                       
[451] "alsfrs.normal_numeric"                                            "alsfrs.slow_or_sloppy_all_words_are_legible_numeric"             
[453] "alsfrs.not_all_words_are_legible_numeric"                         "alsfrs.able_to_grip_pen_numeric"                                 
[455] "alsfrs.unable_to_grip_pen_numeric"                                "alsfrs.have_you_noticed_any_difficulty_walking_numeric"          
[457] "alsfrs.do_you_have_any_difficulty_climbing_stairs_numeric"        "ess.felt_excessively_sleepy_past"                                
[459] "ess.sitting_and_reading"                                          "ess.watching_tv"                                                 
[461] "ess.sitting_inactive_in_a_public_place"                           "ess.break_passenger_car_hour"                                    
[463] "ess.afternoon_rest_circumstances_lying"                           "ess.sitting_and_talking_to_someone"                              
[465] "ess.sitting_quietly_after_a_lunch_without_alcohol"                "ess.car_minutes_stopped_traffic"                                 
[467] "ess.insomnia"                                                     "ess.narcolepsy"                                                  
[469] "ess.rem_sleep_behaviour_disorder"                                 "ess.restless_legs_syndrome"                                      
[471] "ess.sleep_apnea"                                                  "ess.sleep_terrors"                                               
[473] "ess.sleepwalking"                                                 "ess.other"                                                       
[475] "ess.none"                                                         "ess.felt_excessively_sleepy_past_numeric"                        
[477] "ess.sitting_and_reading_numeric"                                  "ess.watching_tv_numeric"                                         
[479] "ess.sitting_inactive_in_a_public_place_numeric"                   "ess.break_passenger_car_hour_numeric"                            
[481] "ess.afternoon_rest_circumstances_lying_numeric"                   "ess.sitting_and_talking_to_someone_numeric"                      
[483] "ess.sitting_quietly_after_a_lunch_without_alcohol_numeric"        "ess.car_minutes_stopped_traffic_numeric"                         
[485] "ess.insomnia_numeric"                                             "ess.narcolepsy_numeric"                                          
[487] "ess.rem_sleep_behaviour_disorder_numeric"                         "ess.restless_legs_syndrome_numeric"                              
[489] "ess.sleep_apnea_numeric"                                          "ess.sleep_terrors_numeric"                                       
[491] "ess.sleepwalking_numeric"                                         "ess.other_numeric"                                               
[493] "ess.none_numeric"                                                 "midas.headaches_suffered_past_month"                             
[495] "midas.days_past_month_miss_work_school"                           "midas.days_past_month_work_school_half_productivity"             
[497] "midas.days_past_month_miss_housework"                             "midas.days_past_month_housework_half_productivity"               
[499] "midas.days_past_month_miss_social_activities"                     "midas.headaches_suffered_past_month_numeric"                     
[501] "midas.days_past_month_miss_work_school_numeric"                   "midas.days_past_month_work_school_half_productivity_numeric"     
[503] "midas.days_past_month_miss_housework_numeric"                     "midas.days_past_month_housework_half_productivity_numeric"       
[505] "midas.days_past_month_miss_social_activities_numeric"             "nmss.ability_taste_food_smell"                                   
[507] "nmss.severe_taste_smell_loss"                                     "nmss.ability_taste_food_smell_numeric"                           
[509] "nmss.severe_taste_smell_loss_numeric"                             "nmsw.drink_problems_choking_difficulty"                          
[511] "nmsw.regularly_night_pass_urine"                                  "nmsw.sex_difficult_finding"                                      
[513] "nmsw.weak_standing_feeling_light"                                 "nmsw.acting_moving_sleep_talking"                                
[515] "nmsw.double_vision"                                               "nmsw.drink_problems_choking_difficulty_numeric"                  
[517] "nmsw.regularly_night_pass_urine_numeric"                          "nmsw.sex_difficult_finding_numeric"                              
[519] "nmsw.weak_standing_feeling_light_numeric"                         "nmsw.acting_moving_sleep_talking_numeric"                        
[521] "nmsw.double_vision_numeric"                                       "tinnitus.ears_noises_minutes_head"                               
[523] "tinnitus.ears_noises_minutes_head_frequency"                      "tinnitus.noises_worry_worst_upset"                               
[525] "tinnitus.ears_noises_minutes_head_numeric"                        "tinnitus.ears_noises_minutes_head_frequency_numeric"             
[527] "tinnitus.noises_worry_worst_upset_numeric"                        "ID.dup"                                                          
[529] "prepandemic_january_31_numeric"                                   "prepandemic_january_31"                                          
[531] "prepandemic_march_1_numeric"                                      "prepandemic_march_1"                                             
[533] "prepandemic_march_23_numeric"                                     "prepandemic_march_23"                                            
[535] "time_diff_coping_lockdown_march_23"                               "time_diff_coping_awareness_march_1"                              
[537] "time_diff_coping_first_case_uk_jan_31"                            "time_diff_signup_lockdown_march_23"                              
[539] "time_diff_signup_awareness_march_1"                               "time_diff_signup_first_case_uk_jan_31"                           
[541] "time_diff_sign_up_coping"                                         "phq.diff_score_retro_prepan"                                     
[543] "phq.diff_score_base_prepan"                                       "phq.diff_score_base_retro"                                       
[545] "gad.diff_score_retro_prepan"                                      "gad.diff_score_base_prepan"                                      
[547] "gad.diff_score_base_retro"                                        "pcl.diff_score_base_prepan"                                      
[549] "ocir.diff_score_retro_prepan"                                     "ocir.diff_score_base_prepan"                                     
[551] "ocir.diff_score_base_retro"                                       "data_group_phq_numeric"                                          
[553] "data_group_phq"                                                   "data_group_gad_numeric"                                          
[555] "data_group_gad"                                                   "data_group_pcl_numeric"                                          
[557] "data_group_pcl"                                                   "data_group_ocir_numeric"                                         
[559] "data_group_ocir"                                                  "respiratory.breath_diarrhoea_fatigue_fever"                      
[561] "respiratory.worst_ill_feel"                                       "respiratory.persistent_cough"                                    
[563] "respiratory.skipping_meals_because_you_felt_unwell"               "respiratory.was_your_sense_of_smell_impaired"                    
[565] "respiratory.was_your_sense_of_taste_impaired"                     "respiratory.feeling_tiredfatigued"                               
[567] "respiratory.throat_swab_coronavirus_nose"                         "respiratory.throat_swab_test_nose"                               
[569] "respiratory.coronavirus_antibody_test"                            "respiratory.results_antibody_test"                               
[571] "respiratory.breath_diarrhoea_fatigue_fever_numeric"               "respiratory.worst_ill_feel_numeric"                              
[573] "respiratory.persistent_cough_numeric"                             "respiratory.skipping_meals_because_you_felt_unwell_numeric"      
[575] "respiratory.was_your_sense_of_smell_impaired_numeric"             "respiratory.was_your_sense_of_taste_impaired_numeric"            
[577] "respiratory.feeling_tiredfatigued_numeric"                        "respiratory.throat_swab_coronavirus_nose_numeric"                
[579] "respiratory.throat_swab_test_nose_numeric"                        "respiratory.coronavirus_antibody_test_numeric"                   
[581] "respiratory.results_antibody_test_numeric"                        "respiratory.anosmia_binary"                                      
[583] "respiratory.pcough_binary"                                        "respiratory.fatigue_binary"                                      
[585] "respiratory.skipped_meals_binary"                                 "Sex_respiratory"                                                 
[587] "Age_respiratory"                                                  "covid_symptoms"                                                  
[589] "covid_probability"                                                "covid_odds_ratio"                                                
[591] "covid_odds_ratio_50_numeric"                                      "covid_odds_ratio_50"                                             
[593] "covid_probable_case_incl_screening_numeric"                       "covid_probable_case_incl_screening"                              
[595] "COPING_ID"                                                        "Birthyear_unc"                                                   
[597] "Birthmonth"                                                       "Birthday"                                                        
[599] "Birthyear_unc_numeric"                                            "Birthmonth_numeric"                                              
[601] "Birthday_numeric"                                                 "Birthyear_numeric"                                               
[603] "Birthdate"                                                        "Cohort"                                                          
[605] "Ethnicity.unc"                                                    "Cohort_numeric"                                                  
[607] "Ethnicity.unc_numeric"                                            "Ethnicity_full"                                                  
[609] "Transgender_uncleaned"                                            "Transgender_uncleaned_numeric"                                   
[611] "employment.what_is_your_highest_level_of_education"               "employment.what_is_your_highest_level_of_education_numeric"      
[613] "respiratory.coronavirus_week_antibody_test"                       "respiratory.coronavirus_week_antibody_test_numeric"              
[615] "covid_test_positive_double_numeric"                               "covid_test_positive_double"                                      
[617] "gad.covid_test_positive_numeric"                                  "covid_test_positive"                                             
[619] "age_category_collapsed_numeric"                                   "age_category_collapsed"                                          
[621] "age_category_collapsed_reordered"                                 "Gender_collapsed_numeric"                                        
[623] "Gender_collapsed"                                                 "Ethnicity_collapsed_numeric"                                     
[625] "Ethnicity_collapsed"                                              "Disorder_hierarchical_psychosis_numeric"                         
[627] "Disorder_hierarchical_psychosis"                                  "Disorder_hierarchical_numeric"                                   
[629] "Disorder_hierarchical"                                            "gad.change_binary_numeric"                                       
[631] "gad.change_binary"                                                "phq.change_binary_numeric"                                       
[633] "phq.change_binary"                                                "ocir.change_binary_numeric"                                      
[635] "ocir.change_binary"                                               "prepandemic_3_months_numeric"                                    
[637] "prepandemic_3_months"                                             "prepandemic_6_months_numeric"                                    
[639] "prepandemic_6_months"                                             "prepandemic_9_months_numeric"                                    
[641] "prepandemic_9_months"                                            

Important vectors with collections of columns

Vectors for sum_score

 [1] "gad.sum_score_prepan"        "pcl.sum_score_prepan"        "phq.sum_score_prepan"        "phq.sum_score_8items_prepan"
 [5] "ocir.sum_score_prepan"       "gad.sum_score_base"          "gad.sum_score_retro"         "pcl.sum_score_base"         
 [9] "pcl.sum_score_retro"         "phq.sum_score_base"          "phq.sum_score_8items_base"   "phq.sum_score_retro"        
[13] "phq.sum_score_8items_retro"  "ocir.sum_score_base"         "ocir.sum_score_retro"       

[1] "phq.sum_score_prepan" "phq.sum_score_base"   "phq.sum_score_retro" 

[1] "gad.sum_score_prepan" "gad.sum_score_base"   "gad.sum_score_retro" 

[1] "ocir.sum_score_prepan" "ocir.sum_score_base"   "ocir.sum_score_retro" 

[1] "pcl.sum_score_prepan" "pcl.sum_score_base"   "pcl.sum_score_retro" 

Vectors for diff_score

 [1] "phq.diff_score_retro_prepan"  "phq.diff_score_base_prepan"   "phq.diff_score_base_retro"    "gad.diff_score_retro_prepan" 
 [5] "gad.diff_score_base_prepan"   "gad.diff_score_base_retro"    "pcl.diff_score_base_prepan"   "ocir.diff_score_retro_prepan"
 [9] "ocir.diff_score_base_prepan"  "ocir.diff_score_base_retro"  

[1] "phq.diff_score_retro_prepan" "phq.diff_score_base_prepan"  "phq.diff_score_base_retro"  

[1] "gad.diff_score_retro_prepan" "gad.diff_score_base_prepan"  "gad.diff_score_base_retro"  

[1] "ocir.diff_score_retro_prepan" "ocir.diff_score_base_prepan"  "ocir.diff_score_base_retro"  

#ANALYSIS PLAN 1. Calculate mean change per disorder group i. Disorders (9) a. Rank the disorders by the size of the change ii. Sex iii. Age (13; RAMP; potentially collapse above 65 depending on group sizes) iv. Ethnicity groups (5; White European, Black (British), Asian, Arab, Mixed)

1. Questionnaire comparions: Calculate mean change per group

gad

Summary scores

Warning in min(., na.rm = T): no non-missing arguments to min; returning Inf

Warning in min(., na.rm = T): no non-missing arguments to min; returning Inf

Warning in max(., na.rm = T): no non-missing arguments to max; returning -Inf

Warning in max(., na.rm = T): no non-missing arguments to max; returning -Inf

# A tibble: 4 x 22
  gad.change_bina… Prepan_mean Base_mean Retro_mean Prepan_sd Base_sd Retro_sd Prepan_median Base_median Retro_median Prepan_min Base_min Retro_min
  <fct>                  <dbl>     <dbl>      <dbl>     <dbl>   <dbl>    <dbl>         <dbl>       <dbl>        <dbl>      <dbl>    <dbl>     <dbl>
1 No change               8.25      2.67     NaN         6.20    4.45    NA                7         0             NA          0        0       Inf
2 Worsening               9.40      8.70       4.64      5.86    5.70     4.28             8         7              4          0        0         0
3 Improvement             9.23      4.77       8.36      5.82    4.83     5.47             8         3              7          0        0         0
4 <NA>                   10         5.7      NaN         8.56    5.99    NA               13         4.5           NA          0        0       Inf
# … with 9 more variables: Prepan_max <dbl>, Base_max <dbl>, Retro_max <dbl>, Prepan_skew <dbl>, Base_skew <dbl>, Retro_skew <dbl>,
#   Prepan_kurtosis <dbl>, Base_kurtosis <dbl>, Retro_kurtosis <dbl>

# A tibble: 12 x 9
   Change      Variable   Min   Max   Mean Median    SD    Skew Kurtosis
   <fct>       <chr>    <dbl> <dbl>  <dbl>  <dbl> <dbl>   <dbl>    <dbl>
 1 No change   Base         0    21   2.67    0    4.45  2.18      4.45 
 2 No change   Prepan       0    21   8.25    7    6.20  0.475    -0.903
 3 No change   Retro      Inf  -Inf NaN      NA   NA    NA        NA    
 4 Worsening   Base         0    21   8.70    7    5.70  0.525    -0.801
 5 Worsening   Prepan       0    21   9.40    8    5.86  0.339    -0.933
 6 Worsening   Retro        0    21   4.64    4    4.28  1.38      1.93 
 7 Improvement Base         0    21   4.77    3    4.83  1.37      1.57 
 8 Improvement Prepan       0    21   9.23    8    5.82  0.330    -0.903
 9 Improvement Retro        0    21   8.36    7    5.47  0.512    -0.587
10 <NA>        Base         0    21   5.7     4.5  5.99  1.04      0.101
11 <NA>        Prepan       0    21  10      13    8.56 -0.0314   -1.99 
12 <NA>        Retro      Inf  -Inf NaN      NA   NA    NA        NA    

Table for knitting gad sum score summary table

Differences scores

Warning in min(., na.rm = T): no non-missing arguments to min; returning Inf

Warning in min(., na.rm = T): no non-missing arguments to min; returning Inf

Warning in min(., na.rm = T): no non-missing arguments to min; returning Inf

Warning in min(., na.rm = T): no non-missing arguments to min; returning Inf

Warning in max(., na.rm = T): no non-missing arguments to max; returning -Inf

Warning in max(., na.rm = T): no non-missing arguments to max; returning -Inf

Warning in max(., na.rm = T): no non-missing arguments to max; returning -Inf

Warning in max(., na.rm = T): no non-missing arguments to max; returning -Inf

# A tibble: 4 x 22
  gad.change_bina… `Base - prepan_… `Base - retro_m… `Prepan - retro… `Base - prepan_… `Base - retro_s… `Prepan - retro… `Base - prepan_…
  <fct>                       <dbl>            <dbl>            <dbl>            <dbl>            <dbl>            <dbl>            <dbl>
1 No change                  -2.05            NaN             NaN                 4.73            NA               NA                  -1
2 Worsening                   0.969             4.06           -3.58              5.31             4.29             5.11                1
3 Improvement                -4.08             -3.59            0.354             5.49             5.32             5.69               -3
4 <NA>                       -1.5             NaN             NaN                 5.43            NA               NA                  -2
# … with 14 more variables: `Base - retro_median` <dbl>, `Prepan - retro_median` <dbl>, `Base - prepan_min` <dbl>, `Base - retro_min` <dbl>, `Prepan
#   - retro_min` <dbl>, `Base - prepan_max` <dbl>, `Base - retro_max` <dbl>, `Prepan - retro_max` <dbl>, `Base - prepan_skew` <dbl>, `Base -
#   retro_skew` <dbl>, `Prepan - retro_skew` <dbl>, `Base - prepan_kurtosis` <dbl>, `Base - retro_kurtosis` <dbl>, `Prepan - retro_kurtosis` <dbl>

# A tibble: 12 x 9
   Change      Variable         Min   Max    Mean Median    SD    Skew Kurtosis
   <fct>       <chr>          <dbl> <dbl>   <dbl>  <dbl> <dbl>   <dbl>    <dbl>
 1 No change   Base - prepan    -21    20  -2.05      -1  4.73 -0.487     1.51 
 2 No change   Base - retro     Inf  -Inf NaN         NA NA    NA        NA    
 3 No change   Prepan - retro   Inf  -Inf NaN         NA NA    NA        NA    
 4 Worsening   Base - prepan    -18    21   0.969      1  5.31 -0.0146    0.601
 5 Worsening   Base - retro     -21    21   4.06       3  4.29  0.770     0.879
 6 Worsening   Prepan - retro   -21    18  -3.58      -3  5.11 -0.289     0.441
 7 Improvement Base - prepan    -21    12  -4.08      -3  5.49 -0.297     0.126
 8 Improvement Base - retro     -21    21  -3.59      -3  5.32  0.475     2.15 
 9 Improvement Prepan - retro   -21    20   0.354      1  5.69 -0.158     0.875
10 <NA>        Base - prepan     -9     6  -1.5       -2  5.43  0.0468   -1.68 
11 <NA>        Base - retro     Inf  -Inf NaN         NA NA    NA        NA    
12 <NA>        Prepan - retro   Inf  -Inf NaN         NA NA    NA        NA    

Table for knitting gad diff summary table

Long format

  ID Sample Gender_collapsed Ethnicity_collapsed Disorder_hierarchical gad.diff_score   Time_point
1  8   GLAD             Male            European  Depressive disorders             NA retro_prepan
2 10   GLAD           Female            European  Depressive disorders             -5 retro_prepan
3 11   GLAD             Male            European  Depressive disorders              7 retro_prepan
4 13   GLAD           Female            European  Depressive disorders              0 retro_prepan
5 15   GLAD           Female            European  Depressive disorders             NA retro_prepan
6 20   GLAD           Female            European  Depressive disorders             15 retro_prepan

  ID Sample Gender_collapsed Ethnicity_collapsed Disorder_hierarchical gad.sum_score Time_point Full_cohort Time_point_dup
1  8   GLAD             Male            European  Depressive disorders            16     prepan Full cohort         prepan
2 10   GLAD           Female            European  Depressive disorders             7     prepan Full cohort         prepan
3 11   GLAD             Male            European  Depressive disorders             2     prepan Full cohort         prepan
4 13   GLAD           Female            European  Depressive disorders             5     prepan Full cohort         prepan
5 15   GLAD           Female            European  Depressive disorders             3     prepan Full cohort         prepan
6 20   GLAD           Female            European  Depressive disorders             6     prepan Full cohort         prepan

Plots diff

makes this into percentages

Warning: Removed 56536 rows containing non-finite values (stat_count).

Warning: Removed 33099 rows containing non-finite values (stat_bin).

Warning: Removed 33099 rows containing non-finite values (stat_bin).

Warning: Removed 23437 rows containing non-finite values (stat_bin).

Plots sum

Warning: Removed 33099 rows containing non-finite values (stat_bin).

Warning: Removed 33099 rows containing non-finite values (stat_bin).

Violin plots

Warning: `fun.y` is deprecated. Use `fun` instead.

Warning: `fun.ymin` is deprecated. Use `fun.min` instead.

Warning: `fun.ymax` is deprecated. Use `fun.max` instead.

Warning: Removed 33099 rows containing non-finite values (stat_ydensity).

Warning: Removed 33099 rows containing non-finite values (stat_summary).

Warning: Removed 33099 rows containing missing values (geom_point).

Warning: `fun.y` is deprecated. Use `fun` instead.

Warning: `fun.ymin` is deprecated. Use `fun.min` instead.

Warning: `fun.ymax` is deprecated. Use `fun.max` instead.

Warning: Removed 33099 rows containing non-finite values (stat_ydensity).

Warning: Removed 33099 rows containing non-finite values (stat_summary).

Warning: Removed 33099 rows containing missing values (geom_point).

Warning: `fun.y` is deprecated. Use `fun` instead.

Warning: `fun.ymin` is deprecated. Use `fun.min` instead.

Warning: `fun.ymax` is deprecated. Use `fun.max` instead.

Warning: Removed 33099 rows containing non-finite values (stat_ydensity).

Warning: Removed 33099 rows containing non-finite values (stat_summary).

Warning: Removed 33099 rows containing missing values (geom_point).

Warning: `fun.y` is deprecated. Use `fun` instead.

Warning: `fun.ymin` is deprecated. Use `fun.min` instead.

Warning: `fun.ymax` is deprecated. Use `fun.max` instead.

Warning: Removed 33099 rows containing non-finite values (stat_ydensity).

Warning: Removed 33099 rows containing non-finite values (stat_summary).

Warning: Removed 33099 rows containing missing values (geom_point).

phq

Summary scores

Warning in min(., na.rm = T): no non-missing arguments to min; returning Inf

Warning in min(., na.rm = T): no non-missing arguments to min; returning Inf

Warning in max(., na.rm = T): no non-missing arguments to max; returning -Inf

Warning in max(., na.rm = T): no non-missing arguments to max; returning -Inf

# A tibble: 4 x 22
  phq.change_bina… Prepan_mean Base_mean Retro_mean Prepan_sd Base_sd Retro_sd Prepan_median Base_median Retro_median Prepan_min Base_min Retro_min
  <fct>                  <dbl>     <dbl>      <dbl>     <dbl>   <dbl>    <dbl>         <dbl>       <dbl>        <dbl>      <dbl>    <dbl>     <dbl>
1 No change               10.5      3.78     NaN         7.21    5.48    NA                9           2           NA          0        0       Inf
2 Worsening               11.7     10.3        5.71      6.74    6.60     5.28            11           9            4          0        0         0
3 Improvement             11.5      6.76       9.81      6.72    5.76     6.78            11           5            9          0        0         0
4 <NA>                    11.4      5.88     NaN         7.58    6.89    NA                9           3           NA          0        0       Inf
# … with 9 more variables: Prepan_max <dbl>, Base_max <dbl>, Retro_max <dbl>, Prepan_skew <dbl>, Base_skew <dbl>, Retro_skew <dbl>,
#   Prepan_kurtosis <dbl>, Base_kurtosis <dbl>, Retro_kurtosis <dbl>

# A tibble: 12 x 9
   Phq.change_binary Variable   Min   Max   Mean Median    SD   Skew Kurtosis
   <fct>             <chr>    <dbl> <dbl>  <dbl>  <dbl> <dbl>  <dbl>    <dbl>
 1 No change         Base         0    27   3.78      2  5.48  2.02     3.82 
 2 No change         Prepan       0    27  10.5       9  7.21  0.474   -0.748
 3 No change         Retro      Inf  -Inf NaN        NA NA    NA       NA    
 4 Worsening         Base         0    27  10.3       9  6.60  0.593   -0.541
 5 Worsening         Prepan       0    27  11.7      11  6.74  0.303   -0.773
 6 Worsening         Retro        0    27   5.71      4  5.28  1.42     1.95 
 7 Improvement       Base         0    27   6.76      5  5.76  1.06     0.707
 8 Improvement       Prepan       0    27  11.5      11  6.72  0.305   -0.812
 9 Improvement       Retro        0    27   9.81      9  6.78  0.596   -0.462
10 <NA>              Base         0    27   5.88      3  6.89  1.15     0.368
11 <NA>              Prepan       0    25  11.4       9  7.58  0.385   -1.14 
12 <NA>              Retro      Inf  -Inf NaN        NA NA    NA       NA    

Table for knitting phq sum score summary table

Differences score

Warning in min(., na.rm = T): no non-missing arguments to min; returning Inf

Warning in min(., na.rm = T): no non-missing arguments to min; returning Inf

Warning in min(., na.rm = T): no non-missing arguments to min; returning Inf

Warning in min(., na.rm = T): no non-missing arguments to min; returning Inf

Warning in max(., na.rm = T): no non-missing arguments to max; returning -Inf

Warning in max(., na.rm = T): no non-missing arguments to max; returning -Inf

Warning in max(., na.rm = T): no non-missing arguments to max; returning -Inf

Warning in max(., na.rm = T): no non-missing arguments to max; returning -Inf

# A tibble: 4 x 22
  phq.change_bina… `Base - prepan_… `Base - retro_m… `Prepan - retro… `Base - prepan_… `Base - retro_s… `Prepan - retro… `Base - prepan_…
  <fct>                       <dbl>            <dbl>            <dbl>            <dbl>            <dbl>            <dbl>            <dbl>
1 No change                  -1.98            NaN             NaN                 5.13            NA               NA                  -1
2 Worsening                   0.800             4.60           -4.25              5.55             4.54             5.50                1
3 Improvement                -4.03             -3.05           -0.121             5.60             5.58             6.13               -3
4 <NA>                       -0.531           NaN             NaN                 6.11            NA               NA                   0
# … with 14 more variables: `Base - retro_median` <dbl>, `Prepan - retro_median` <dbl>, `Base - prepan_min` <dbl>, `Base - retro_min` <dbl>, `Prepan
#   - retro_min` <dbl>, `Base - prepan_max` <dbl>, `Base - retro_max` <dbl>, `Prepan - retro_max` <dbl>, `Base - prepan_skew` <dbl>, `Base -
#   retro_skew` <dbl>, `Prepan - retro_skew` <dbl>, `Base - prepan_kurtosis` <dbl>, `Base - retro_kurtosis` <dbl>, `Prepan - retro_kurtosis` <dbl>

# A tibble: 12 x 9
   Phq.change_binary Variable         Min   Max    Mean Median    SD    Skew Kurtosis
   <fct>             <chr>          <dbl> <dbl>   <dbl>  <dbl> <dbl>   <dbl>    <dbl>
 1 No change         Base - prepan    -25    23  -1.98      -1  5.13 -0.317     1.69 
 2 No change         Base - retro     Inf  -Inf NaN         NA NA    NA        NA    
 3 No change         Prepan - retro   Inf  -Inf NaN         NA NA    NA        NA    
 4 Worsening         Base - prepan    -23    24   0.800      1  5.55 -0.0686    0.878
 5 Worsening         Base - retro     -18    27   4.60       4  4.54  0.927     1.34 
 6 Worsening         Prepan - retro   -27    21  -4.25      -4  5.50 -0.377     0.604
 7 Improvement       Base - prepan    -22    14  -4.03      -3  5.60 -0.412     0.231
 8 Improvement       Base - retro     -24    23  -3.05      -3  5.58  0.0340    1.85 
 9 Improvement       Prepan - retro   -27    22  -0.121      0  6.13 -0.0254    0.730
10 <NA>              Base - prepan    -21     8  -0.531      0  6.11 -1.10      1.88 
11 <NA>              Base - retro     Inf  -Inf NaN         NA NA    NA        NA    
12 <NA>              Prepan - retro   Inf  -Inf NaN         NA NA    NA        NA    

Table for knitting phq diff summary table

Long format

  ID Sample Gender_collapsed Ethnicity_collapsed Disorder_hierarchical phq.diff_score   Time_point
1  8   GLAD             Male            European  Depressive disorders              9 retro_prepan
2 10   GLAD           Female            European  Depressive disorders             -3 retro_prepan
3 11   GLAD             Male            European  Depressive disorders              5 retro_prepan
4 13   GLAD           Female            European  Depressive disorders             -1 retro_prepan
5 15   GLAD           Female            European  Depressive disorders             NA retro_prepan
6 20   GLAD           Female            European  Depressive disorders              7 retro_prepan

  ID Sample Gender_collapsed Ethnicity_collapsed Disorder_hierarchical phq.sum_score Time_point
1  8   GLAD             Male            European  Depressive disorders            14     prepan
2 10   GLAD           Female            European  Depressive disorders             4     prepan
3 11   GLAD             Male            European  Depressive disorders             2     prepan
4 13   GLAD           Female            European  Depressive disorders             8     prepan
5 15   GLAD           Female            European  Depressive disorders             5     prepan
6 20   GLAD           Female            European  Depressive disorders            13     prepan

Plots

makes this into percentages

Warning: Removed 55288 rows containing non-finite values (stat_count).

Warning: Removed 32162 rows containing non-finite values (stat_bin).

Warning: Removed 23126 rows containing non-finite values (stat_bin).

Warning: Removed 32162 rows containing non-finite values (stat_bin).

Plots sum

Warning: Removed 32162 rows containing non-finite values (stat_bin).

Warning: Removed 32162 rows containing non-finite values (stat_bin).

Violin plots

Warning: `fun.y` is deprecated. Use `fun` instead.

Warning: `fun.ymin` is deprecated. Use `fun.min` instead.

Warning: `fun.ymax` is deprecated. Use `fun.max` instead.

Warning: Removed 32162 rows containing non-finite values (stat_ydensity).

Warning: Removed 32162 rows containing non-finite values (stat_summary).

Warning: Removed 32162 rows containing missing values (geom_point).

Warning: `fun.y` is deprecated. Use `fun` instead.

Warning: `fun.ymin` is deprecated. Use `fun.min` instead.

Warning: `fun.ymax` is deprecated. Use `fun.max` instead.

Warning: Removed 32162 rows containing non-finite values (stat_ydensity).

Warning: Removed 32162 rows containing non-finite values (stat_summary).

Warning: Removed 32162 rows containing missing values (geom_point).

Warning: `fun.y` is deprecated. Use `fun` instead.

Warning: `fun.ymin` is deprecated. Use `fun.min` instead.

Warning: `fun.ymax` is deprecated. Use `fun.max` instead.

Warning: Removed 32162 rows containing non-finite values (stat_ydensity).

Warning: Removed 32162 rows containing non-finite values (stat_summary).

Warning: Removed 32162 rows containing missing values (geom_point).

Warning: `fun.y` is deprecated. Use `fun` instead.

Warning: `fun.ymin` is deprecated. Use `fun.min` instead.

Warning: `fun.ymax` is deprecated. Use `fun.max` instead.

Warning: Removed 32162 rows containing non-finite values (stat_ydensity).

Warning: Removed 32162 rows containing non-finite values (stat_summary).

Warning: Removed 32162 rows containing missing values (geom_point).

ocir

Summary scores

Warning in min(., na.rm = T): no non-missing arguments to min; returning Inf

Warning in max(., na.rm = T): no non-missing arguments to max; returning -Inf

# A tibble: 4 x 22
  ocir.change_bin… Prepan_mean Base_mean Retro_mean Prepan_sd Base_sd Retro_sd Prepan_median Base_median Retro_median Prepan_min Base_min Retro_min
  <fct>                  <dbl>     <dbl>      <dbl>     <dbl>   <dbl>    <dbl>         <dbl>       <dbl>        <dbl>      <dbl>    <dbl>     <dbl>
1 No change               25.4      14.3       17.4      13.0    12.3     12.7          24            11           15          0        0         0
2 Worsening               27.5      17.8       17.0      14.1    12.9     12.8          26            15           14          1        0         0
3 Improvement             27.5      21.8       18.1      14.5    14.1     11.5          23.5          18           16         10        0         0
4 <NA>                    16.4     NaN         17        10.7    NA       NA            14            NA           17          0      Inf        17
# … with 9 more variables: Prepan_max <dbl>, Base_max <dbl>, Retro_max <dbl>, Prepan_skew <dbl>, Base_skew <dbl>, Retro_skew <dbl>,
#   Prepan_kurtosis <dbl>, Base_kurtosis <dbl>, Retro_kurtosis <dbl>

# A tibble: 3 x 8
  Variable   Min   Max  Mean Median    SD  Skew Kurtosis
  <chr>    <dbl> <dbl> <dbl>  <dbl> <dbl> <dbl>    <dbl>
1 Base         0    72  17.8     15  12.9 1.13    1.17  
2 Prepan       1    71  27.5     26  14.1 0.627  -0.0287
3 Retro        0    71  17.0     14  12.8 1.10    0.965 

Table for knitting ocir sum score summary table

Differences scores

  Base - prepan_mean Base - retro_mean Prepan - retro_mean Base - prepan_sd Base - retro_sd Prepan - retro_sd Base - prepan_median
1          -8.043956         0.9100517             -3.5058         10.74915        9.658776          9.301472                   -7
  Base - retro_median Prepan - retro_median Base - prepan_min Base - retro_min Prepan - retro_min Base - prepan_max Base - retro_max
1                   0                    -3               -40              -44                -30                34               72
  Prepan - retro_max Base - prepan_skew Base - retro_skew Prepan - retro_skew Base - prepan_kurtosis Base - retro_kurtosis Prepan - retro_kurtosis
1                 27         -0.1101172         0.2846682         -0.06193155              0.6746143              2.470417               0.3779813

        Variable Min Max       Mean Median        SD        Skew  Kurtosis
1  Base - prepan -40  34 -8.0439560     -7 10.749154 -0.11011724 0.6746143
2   Base - retro -44  72  0.9100517      0  9.658776  0.28466820 2.4704171
3 Prepan - retro -30  27 -3.5058005     -3  9.301472 -0.06193155 0.3779813

Table for knitting ocir diff summary table

Long format

  ID Sample Gender_collapsed Ethnicity_collapsed Disorder_hierarchical ocir.diff_score   Time_point
1  8   GLAD             Male            European  Depressive disorders              NA retro_prepan
2 10   GLAD           Female            European  Depressive disorders              NA retro_prepan
3 11   GLAD             Male            European  Depressive disorders              NA retro_prepan
4 13   GLAD           Female            European  Depressive disorders              NA retro_prepan
5 15   GLAD           Female            European  Depressive disorders              NA retro_prepan
6 20   GLAD           Female            European  Depressive disorders              NA retro_prepan

  ID Sample Gender_collapsed Ethnicity_collapsed Disorder_hierarchical ocir.sum_score Time_point
1  8   GLAD             Male            European  Depressive disorders             NA     prepan
2 10   GLAD           Female            European  Depressive disorders             NA     prepan
3 11   GLAD             Male            European  Depressive disorders             NA     prepan
4 13   GLAD           Female            European  Depressive disorders             NA     prepan
5 15   GLAD           Female            European  Depressive disorders             NA     prepan
6 20   GLAD           Female            European  Depressive disorders             NA     prepan

Plots

makes this into percentages

Warning: Removed 89888 rows containing non-finite values (stat_count).

Warning: Removed 58866 rows containing non-finite values (stat_bin).

Warning: Removed 31022 rows containing non-finite values (stat_bin).

Warning: Removed 58866 rows containing non-finite values (stat_bin).

Plots sum

Warning: Removed 81548 rows containing non-finite values (stat_bin).

Warning: Removed 81548 rows containing non-finite values (stat_bin).

Violin plots

Warning: `fun.y` is deprecated. Use `fun` instead.

Warning: `fun.ymin` is deprecated. Use `fun.min` instead.

Warning: `fun.ymax` is deprecated. Use `fun.max` instead.

Warning: Removed 81548 rows containing non-finite values (stat_ydensity).

Warning: Removed 81548 rows containing non-finite values (stat_summary).

Warning: Removed 81548 rows containing missing values (geom_point).

Warning: `fun.y` is deprecated. Use `fun` instead.

Warning: `fun.ymin` is deprecated. Use `fun.min` instead.

Warning: `fun.ymax` is deprecated. Use `fun.max` instead.

Warning: Removed 81548 rows containing non-finite values (stat_ydensity).

Warning: Removed 81548 rows containing non-finite values (stat_summary).

Warning: Removed 81548 rows containing missing values (geom_point).

Warning: `fun.y` is deprecated. Use `fun` instead.

Warning: `fun.ymin` is deprecated. Use `fun.min` instead.

Warning: `fun.ymax` is deprecated. Use `fun.max` instead.

Warning: Removed 81548 rows containing non-finite values (stat_ydensity).

Warning: Removed 81548 rows containing non-finite values (stat_summary).

Warning: Removed 81548 rows containing missing values (geom_point).

Warning: `fun.y` is deprecated. Use `fun` instead.

Warning: `fun.ymin` is deprecated. Use `fun.min` instead.

Warning: `fun.ymax` is deprecated. Use `fun.max` instead.

Warning: Removed 81548 rows containing non-finite values (stat_ydensity).

Warning: Removed 81548 rows containing non-finite values (stat_summary).

Warning: Removed 81548 rows containing missing values (geom_point).

++CH: Rank the disorders by the size of the change

PCL: Calculate the total GLAD, COPING, RAMP baseline

PCL

  Base - prepan_mean Base - prepan_sd Base - prepan_median Base - prepan_min Base - prepan_max Base - prepan_skew Base - prepan_kurtosis
1          0.2068164         5.666463                    0               -23                24        -0.07432746              0.6542706

       Variable Min Max      Mean Median       SD        Skew  Kurtosis
1 Base - prepan -23  24 0.2068164      0 5.666463 -0.07432746 0.6542706

Table for knitting pcl sum score summary table

  Base - prepan_mean Base - prepan_sd Base - prepan_median Base - prepan_min Base - prepan_max Base - prepan_skew Base - prepan_kurtosis
1          0.2068164         5.666463                    0               -23                24        -0.07432746              0.6542706

  Variable Min Max     Mean Median       SD      Skew    Kurtosis
1     Base   6  30 12.33784     11 6.051631 0.9357645 -0.01527039
2   Prepan   6  30 15.63577     15 6.120978 0.3799305 -0.73869515

Table for knitting pcl diff summary table

  ID Sample Gender_collapsed Ethnicity_collapsed Disorder_hierarchical pcl.sum_score Time_point
1  8   GLAD             Male            European  Depressive disorders            18     prepan
2 10   GLAD           Female            European  Depressive disorders            22     prepan
3 11   GLAD             Male            European  Depressive disorders            11     prepan
4 13   GLAD           Female            European  Depressive disorders            15     prepan
5 15   GLAD           Female            European  Depressive disorders            16     prepan
6 20   GLAD           Female            European  Depressive disorders            16     prepan

For each item present percentages of how say that this is due to the pandemic

# A tibble: 18 x 5
# Groups:   PCL Item [6]
   `PCL Item`                                                    Binary_score     n  prop n_prop       
   <chr>                                                         <fct>        <int> <dbl> <chr>        
 1 pcl.cut_people_feeling_distant_retro_numeric                  0             5970  18.7 5970 (18.7%) 
 2 pcl.cut_people_feeling_distant_retro_numeric                  1            15402  48.3 15402 (48.3%)
 3 pcl.cut_people_feeling_distant_retro_numeric                  Missing      10512  33   10512 (33%)  
 4 pcl.difficulty_concentrating_retro_numeric                    0             9916  31.1 9916 (31.1%) 
 5 pcl.difficulty_concentrating_retro_numeric                    1            10289  32.3 10289 (32.3%)
 6 pcl.difficulty_concentrating_retro_numeric                    Missing      11679  36.6 11679 (36.6%)
 7 pcl.feeling_irritable_or_having_angry_outbursts_retro_numeric 0             8756  27.5 8756 (27.5%) 
 8 pcl.feeling_irritable_or_having_angry_outbursts_retro_numeric 1             9870  31   9870 (31%)   
 9 pcl.feeling_irritable_or_having_angry_outbursts_retro_numeric Missing      13258  41.6 13258 (41.6%)
10 pcl.stressful_experience_repeated_images_retro_numeric        0            13112  41.1 13112 (41.1%)
11 pcl.stressful_experience_repeated_images_retro_numeric        1             2902   9.1 2902 (9.1%)  
12 pcl.stressful_experience_repeated_images_retro_numeric        Missing      15870  49.8 15870 (49.8%)
13 pcl.stressful_experience_upset_reminded_retro_numeric         0            13848  43.4 13848 (43.4%)
14 pcl.stressful_experience_upset_reminded_retro_numeric         1             3742  11.7 3742 (11.7%) 
15 pcl.stressful_experience_upset_reminded_retro_numeric         Missing      14294  44.8 14294 (44.8%)
16 pcl.stressful_situation_avoiding_activities_retro_numeric     0             9158  28.7 9158 (28.7%) 
17 pcl.stressful_situation_avoiding_activities_retro_numeric     1             2465   7.7 2465 (7.7%)  
18 pcl.stressful_situation_avoiding_activities_retro_numeric     Missing      20261  63.5 20261 (63.5%)

Plots sum

Warning: Removed 43815 rows containing non-finite values (stat_bin).

Warning: Removed 43815 rows containing non-finite values (stat_bin).

Create subsets of data for analyses

Filter data frame to create 2 additional datasets. Additional exclusion criteria (beyond those listed above) are specified below:

1. Retrospective - baseline change anlyses i. Exclude individuals missing retrospective and/or baseline data
2. Prepandemic - baseline change analyses i. Exclude individuals missing pandemic and/or baseline data ii. Exclude individuals with prepandemic data collected after 1st March.

Retro-baseline analyses

[1] 23539   641

Prepan-baseline analyses Exclude individuals with NA for MHD, age, gender; exclude if prepan data collected after 1 March

[1] 11230   641

Prepan-baseline analyses Exclude individuals with NA for MHD, age, gender; exclude if prepan data collected after 1 March (to ensure comparability with prepan analyses)

[1] 9745  641

age² sensitivity Exclude individuals with NA for continuous Age variable

[1] 24816   641

[1] 17092   641

[1] 11230   641

Sensitivity analysis for individuals who signed up to GLAD (prepan) either 3, 6 or 9 months before the pendemic began.

1. Regression models (Regress covariates out of the change? Which?) i. Model 0: Disorders ii. Model 1: Model 0 + Sex + Age + Ethnicity + Time registration and baseline

2. Re-calculate change

3. Re-rank the change again

4. PCL: i. Calculate the total GLAD, COPING, RAMP baseline ii. For each item: a. present percentages of how say that this is due to the pandemic b. Calculate a variable PCL pandemic score c. If one item is highly endorsed due to pandemic -> exclude in sensitivity analysis

2. Regression

Regress covariates out of the change? Which?

Model naming conventions: 0 - unadjusted 1 - adjusted 2 - adjusted with age² a - hierarchical disorders b - overlapping disorder categories

++MD: We may want to exclude disorders from Model b that aren't including in the hierarchy (e.g. ADHD, ASD)

Functions

GAD baseline, not controlling for prior symptoms

Model 0a: Disorder_hierarchical, Unadjusted

#SB Hashing out this chunk of code as it is not included in our final analysis and is currently breaking the code.

# A tibble: 7 x 6
  `Indepent variable`  Estimate `95% CI low` `95% CI up` SE    `p value`
  <fct>                <chr>    <chr>        <chr>       <chr> <chr>    
1 Intercept            2.48     2.38         2.57        0.05  0.00e+00 
2 Anxiety disorders    4.13     3.88         4.39        0.13  5.84e-217
3 Eating disorders     3.76     2.94         4.58        0.42  3.03e-19 
4 OCDs                 6.52     5.85         7.19        0.34  1.26e-80 
5 PTSD                 4.97     4.31         5.64        0.34  3.36e-48 
6 Depressive disorders 5.76     5.63         5.88        0.06  0.00e+00 
7 Bipolar disorders    7.40     7.08         7.72        0.16  0.00e+00 

# A tibble: 1 x 6
  n     R2    `adjust. R2` F       `p value`    df
  <chr> <chr> <chr>        <chr>   <chr>     <dbl>
1 31331 0.22  0.22         1503.62 0.00e+00      6

Model 0b: Overlapping disorders, Unadjusted

# A tibble: 13 x 6
   `Indepent variable`             Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                       2.67     " 2.59"      2.76        0.04  0.00e+00 
 2 Depressive and anxiety disorder 5.63     " 5.50"      5.77        0.07  0.00e+00 
 3 Only depressive disorder        2.08     " 1.90"      2.25        0.09  4.32e-117
 4 Only anxiety disorder           4.03     " 3.81"      4.26        0.11  1.97e-264
 5 Eating disorders                1.46     " 1.22"      1.69        0.12  1.10e-33 
 6 OCDs                            1.94     " 1.71"      2.17        0.12  8.53e-63 
 7 Psychotic and bipolar disorder  0.38     "-0.26"      1.02        0.33  0.24     
 8 Only psychotic disorder         0.35     "-0.20"      0.90        0.28  0.21     
 9 Only bipolar disorder           0.46     " 0.11"      0.80        0.18  0.009    
10 PTSD                            1.56     " 1.35"      1.77        0.11  2.35e-48 
11 ASD                             2.18     " 1.77"      2.59        0.21  1.80e-25 
12 ADHD                            1.17     " 0.63"      1.71        0.28  2.36e-05 
13 Personality disorder            2.19     " 1.89"      2.49        0.15  5.79e-46 

# A tibble: 1 x 6
  n     R2    `adjust. R2` F       `p value`    df
  <chr> <chr> <chr>        <chr>   <chr>     <dbl>
1 31331 0.22  0.22         1503.62 0.00e+00      6

Model 1a: Disorder hierarchical, Adjusted

Model 0a + Gender_collapsed + age_category_collapsed_reordered + Ethnicity_collapsed + Time registration and baseline

# A tibble: 22 x 6
   `Indepent variable`                   Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                                 <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                             " 5.67"  " 5.38"      " 5.96"     0.15  9.14e-313
 2 Anxiety disorders                     " 2.96"  " 2.70"      " 3.22"     0.13  1.10e-109
 3 Eating disorders                      " 2.10"  " 1.30"      " 2.91"     0.41  3.42e-07 
 4 OCDs                                  " 4.01"  " 3.35"      " 4.67"     0.34  1.26e-32 
 5 PTSD                                  " 3.96"  " 3.30"      " 4.63"     0.34  2.24e-31 
 6 Depressive disorders                  " 4.53"  " 4.39"      " 4.66"     0.07  0.00e+00 
 7 Bipolar disorders                     " 6.08"  " 5.76"      " 6.39"     0.16  1.26e-299
 8 Female                                " 1.00"  " 0.86"      " 1.14"     0.07  3.17e-44 
 9 Non-binary/Self-defined               " 1.60"  " 1.04"      " 2.16"     0.29  2.64e-08 
10 age_category_collapsed_reordered16-18 " 3.23"  " 2.85"      " 3.60"     0.19  1.44e-63 
# … with 12 more rows

# A tibble: 1 x 6
  n     R2    `adjust. R2` F       `p value`    df
  <chr> <chr> <chr>        <chr>   <chr>     <dbl>
1 31331 0.22  0.22         1503.62 0.00e+00      6

Model 1b: Overlapping disorders, Adjusted

Model 0a + Gender_collapsed + age_category_collapsed_reordered + Ethnicity_collapsed + Time registration and baseline

# A tibble: 29 x 6
   `Indepent variable`             Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                       " 5.57"  " 5.28"      " 5.85"     0.14  2.40e-317
 2 Depressive and anxiety disorder " 4.60"  " 4.46"      " 4.75"     0.07  0.00e+00 
 3 Only depressive disorder        " 1.81"  " 1.63"      " 1.99"     0.09  2.60e-83 
 4 Only anxiety disorder           " 3.05"  " 2.82"      " 3.28"     0.12  6.22e-145
 5 Eating disorders                " 0.90"  " 0.67"      " 1.13"     0.12  3.61e-14 
 6 OCDs                            " 1.45"  " 1.23"      " 1.67"     0.11  5.23e-37 
 7 Psychotic and bipolar disorder  " 0.74"  " 0.11"      " 1.37"     0.32  0.02     
 8 Only psychotic disorder         " 0.51"  "-0.03"      " 1.05"     0.27  0.06     
 9 Only bipolar disorder           " 0.65"  " 0.31"      " 0.99"     0.17  1.79e-04 
10 PTSD                            " 1.58"  " 1.38"      " 1.79"     0.11  1.05e-50 
# … with 19 more rows

# A tibble: 1 x 6
  n     R2    `adjust. R2` F       `p value`    df
  <chr> <chr> <chr>        <chr>   <chr>     <dbl>
1 31331 0.22  0.22         1503.62 0.00e+00      6

Model 2a: Disorder hierarchical, Adjusted with Age²

# A tibble: 16 x 6
   `Indepent variable`                   Estimate  `95% CI low` `95% CI up` SE    `p value`
   <fct>                                 <chr>     <chr>        <chr>       <chr> <chr>    
 1 Intercept                             "   1.78" "   1.48"    "   2.09"   0.16  4.29e-30 
 2 Anxiety disorders                     "   2.79" "   2.49"    "   3.10"   0.15  3.88e-72 
 3 Eating disorders                      "   2.27" "   1.33"    "   3.20"   0.48  2.04e-06 
 4 OCDs                                  "   4.46" "   3.65"    "   5.26"   0.41  3.29e-27 
 5 PTSD                                  "   4.24" "   3.37"    "   5.12"   0.45  2.15e-21 
 6 Depressive disorders                  "   4.80" "   4.64"    "   4.95"   0.08  0.00e+00 
 7 Bipolar disorders                     "   6.26" "   5.92"    "   6.60"   0.17  6.20e-274
 8 Female                                "   0.71" "   0.56"    "   0.86"   0.08  2.18e-20 
 9 Non-binary/Self-defined               "   1.44" "   0.84"    "   2.05"   0.31  2.84e-06 
10 Age                                   "-248.64" "-259.74"    "-237.54"   5.66  0.00e+00 
11 Age^2                                 "  29.43" "  19.41"    "  39.45"   5.11  8.72e-09 
12 Mixed or multiple ethnic origins      "   0.03" "  -0.46"    "   0.53"   0.25  0.90     
13 Asian or Asian British                "  -0.06" "  -0.71"    "   0.60"   0.33  0.87     
14 African or African British            "   0.02" "  -1.01"    "   1.06"   0.53  0.96     
15 Other ethnicity                       "   1.44" "   0.51"    "   2.38"   0.48  0.002    
16 Time difference: Lockdown to baseline "   0.04" "   0.00"    "   0.08"   0.02  0.03     

# A tibble: 1 x 6
  n     R2    `adjust. R2` F       `p value`    df
  <chr> <chr> <chr>        <chr>   <chr>     <dbl>
1 24462 0.27  0.27         1541.15 0.00e+00      6

Model 2b: Disorder hierarchical, Adjusted with Age²

# A tibble: 23 x 6
   `Indepent variable`             Estimate  `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>     <chr>        <chr>       <chr> <chr>    
 1 Intercept                       "   1.86" "   1.56"    "   2.15"   0.15  1.11e-34 
 2 Depressive and anxiety disorder "   4.91" "   4.75"    "   5.07"   0.08  0.00e+00 
 3 Only depressive disorder        "   1.98" "   1.78"    "   2.19"   0.10  8.03e-79 
 4 Only anxiety disorder           "   3.02" "   2.75"    "   3.29"   0.14  1.46e-104
 5 Eating disorders                "   1.09" "   0.84"    "   1.34"   0.13  2.20e-17 
 6 OCDs                            "   1.38" "   1.13"    "   1.62"   0.12  8.44e-29 
 7 Psychotic and bipolar disorder  "   0.44" "  -0.27"    "   1.15"   0.36  0.23     
 8 Only psychotic disorder         "   0.48" "  -0.10"    "   1.06"   0.30  0.10     
 9 Only bipolar disorder           "   0.64" "   0.29"    "   0.99"   0.18  3.50e-04 
10 PTSD                            "   1.69" "   1.47"    "   1.92"   0.11  3.91e-49 
# … with 13 more rows

# A tibble: 1 x 6
  n     R2    `adjust. R2` F       `p value`    df
  <chr> <chr> <chr>        <chr>   <chr>     <dbl>
1 24462 0.27  0.27         1541.15 0.00e+00      6

GAD prepan

Regression: Change scores

Model 0a: Disorder_hierarchical, Unadjusted

# A tibble: 7 x 6
  `Indepent variable`  Estimate `95% CI low` `95% CI up` SE    `p value`
  <fct>                <chr>    <chr>        <chr>       <chr> <chr>    
1 Intercept            "-1.16"  -2.09        "-0.22"     0.48  0.02     
2 Anxiety disorders    " 0.19"  -0.84        " 1.22"     0.53  0.72     
3 Eating disorders     "-0.57"  -2.60        " 1.47"     1.04  0.59     
4 OCDs                 "-0.14"  -1.55        " 1.26"     0.72  0.84     
5 PTSD                 "-0.06"  -1.68        " 1.57"     0.83  0.95     
6 Depressive disorders " 0.85"  -0.09        " 1.79"     0.48  0.08     
7 Bipolar disorders    " 0.80"  -0.21        " 1.81"     0.52  0.12     

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 11144 0.00  0.00         2.98  0.007         6

Model 0b: Overlapping disorders, Unadjusted

# A tibble: 13 x 6
   `Indepent variable`             Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                       "-0.83"  "-1.51"      "-0.16"     0.34  0.02     
 2 Depressive and anxiety disorder " 0.45"  "-0.24"      " 1.13"     0.35  0.20     
 3 Only depressive disorder        " 0.98"  " 0.24"      " 1.72"     0.38  0.009    
 4 Only anxiety disorder           "-0.24"  "-1.01"      " 0.53"     0.40  0.54     
 5 Eating disorders                " 0.18"  "-0.15"      " 0.51"     0.17  0.29     
 6 OCDs                            "-0.26"  "-0.57"      " 0.05"     0.16  0.11     
 7 Psychotic and bipolar disorder  "-0.51"  "-1.42"      " 0.39"     0.46  0.27     
 8 Only psychotic disorder         " 0.76"  " 0.00"      " 1.52"     0.39  0.05     
 9 Only bipolar disorder           " 0.20"  "-0.26"      " 0.65"     0.23  0.40     
10 PTSD                            " 0.14"  "-0.16"      " 0.43"     0.15  0.36     
11 ASD                             " 0.16"  "-0.44"      " 0.76"     0.31  0.61     
12 ADHD                            "-0.74"  "-1.56"      " 0.09"     0.42  0.08     
13 Personality disorder            "-0.16"  "-0.57"      " 0.24"     0.21  0.43     

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 11144 0.00  0.00         2.98  0.007         6

Model 1a: Disorder hierarchical, Adjusted

Model 0a + Gender_collapsed + age_category_collapsed_reordered + Ethnicity_collapsed + Time registration and baseline

# A tibble: 22 x 6
   `Indepent variable`                   Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                                 <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                             "-1.46"  "-2.56"      "-0.37"     0.56  0.009    
 2 Anxiety disorders                     " 0.13"  "-0.91"      " 1.16"     0.53  0.81     
 3 Eating disorders                      "-0.72"  "-2.75"      " 1.32"     1.04  0.49     
 4 OCDs                                  "-0.24"  "-1.65"      " 1.16"     0.72  0.74     
 5 PTSD                                  "-0.21"  "-1.85"      " 1.43"     0.84  0.80     
 6 Depressive disorders                  " 0.75"  "-0.19"      " 1.69"     0.48  0.12     
 7 Bipolar disorders                     " 0.76"  "-0.26"      " 1.77"     0.52  0.14     
 8 Female                                " 0.49"  " 0.21"      " 0.77"     0.14  5.29e-04 
 9 Non-binary/Self-defined               " 0.42"  "-0.37"      " 1.22"     0.40  0.29     
10 age_category_collapsed_reordered16-18 "-0.50"  "-1.16"      " 0.17"     0.34  0.14     
# … with 12 more rows

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 11144 0.00  0.00         2.98  0.007         6

Model 1b: Overlapping disorders, Adjusted

Model 1b + Gender_collapsed + age_category_collapsed_reordered + Ethnicity_collapsed + Time registration and baseline

# A tibble: 29 x 6
   `Indepent variable`             Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                       "-1.18"  "-2.07"      "-0.29"     0.46  0.010    
 2 Depressive and anxiety disorder " 0.38"  "-0.31"      " 1.08"     0.35  0.27     
 3 Only depressive disorder        " 0.92"  " 0.18"      " 1.67"     0.38  0.02     
 4 Only anxiety disorder           "-0.29"  "-1.07"      " 0.49"     0.40  0.47     
 5 Eating disorders                " 0.14"  "-0.20"      " 0.47"     0.17  0.42     
 6 OCDs                            "-0.26"  "-0.58"      " 0.05"     0.16  0.10     
 7 Psychotic and bipolar disorder  "-0.40"  "-1.31"      " 0.52"     0.47  0.39     
 8 Only psychotic disorder         " 0.79"  " 0.02"      " 1.55"     0.39  0.04     
 9 Only bipolar disorder           " 0.29"  "-0.17"      " 0.74"     0.23  0.22     
10 PTSD                            " 0.11"  "-0.18"      " 0.41"     0.15  0.46     
# … with 19 more rows

Model 2a: Disorder hierarchical, Adjusted with Age²

Model 0a + Gender_collapsed + age²⁺ Ethnicity_collapsed + Time registration and baseline

# A tibble: 16 x 6
   `Indepent variable`                   Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                                 <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                             "-1.44"  " -2.51"     "-0.36"     0.55  0.009    
 2 Anxiety disorders                     " 0.13"  " -0.90"     " 1.17"     0.53  0.80     
 3 Eating disorders                      "-0.70"  " -2.74"     " 1.33"     1.04  0.50     
 4 OCDs                                  "-0.23"  " -1.64"     " 1.17"     0.72  0.75     
 5 PTSD                                  "-0.19"  " -1.83"     " 1.45"     0.84  0.82     
 6 Depressive disorders                  " 0.76"  " -0.18"     " 1.71"     0.48  0.11     
 7 Bipolar disorders                     " 0.76"  " -0.26"     " 1.78"     0.52  0.14     
 8 Female                                " 0.49"  "  0.21"     " 0.77"     0.14  5.24e-04 
 9 Non-binary/Self-defined               " 0.43"  " -0.37"     " 1.22"     0.40  0.29     
10 Age                                   " 0.90"  "-10.29"     "12.09"     5.71  0.87     
11 Age^2                                 "-9.06"  "-20.02"     " 1.89"     5.59  0.10     
12 Mixed or multiple ethnic origins      " 0.30"  " -0.42"     " 1.02"     0.37  0.42     
13 Asian or Asian British                "-0.25"  " -1.42"     " 0.92"     0.60  0.67     
14 African or African British            "-1.26"  " -3.02"     " 0.50"     0.90  0.16     
15 Other ethnicity                       " 0.04"  " -1.08"     " 1.16"     0.57  0.95     
16 Time difference: Lockdown to baseline "-0.01"  " -0.07"     " 0.06"     0.03  0.85     

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 11144 0.00  0.00         2.98  0.007         6

Model 2b: Overlapping disorders, Adjusted with Age²

Model 0a + Gender_collapsed + age² + Ethnicity_collapsed + Time registration and baseline

# A tibble: 23 x 6
   `Indepent variable`             Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                       "-1.19"  " -2.06"     "-0.33"     0.44  0.007    
 2 Depressive and anxiety disorder " 0.41"  " -0.28"     " 1.10"     0.35  0.25     
 3 Only depressive disorder        " 0.95"  "  0.21"     " 1.70"     0.38  0.01     
 4 Only anxiety disorder           "-0.27"  " -1.05"     " 0.51"     0.40  0.50     
 5 Eating disorders                " 0.13"  " -0.20"     " 0.47"     0.17  0.44     
 6 OCDs                            "-0.26"  " -0.58"     " 0.05"     0.16  0.10     
 7 Psychotic and bipolar disorder  "-0.41"  " -1.32"     " 0.50"     0.47  0.38     
 8 Only psychotic disorder         " 0.79"  "  0.03"     " 1.55"     0.39  0.04     
 9 Only bipolar disorder           " 0.28"  " -0.17"     " 0.74"     0.23  0.23     
10 PTSD                            " 0.10"  " -0.19"     " 0.40"     0.15  0.50     
# … with 13 more rows

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 11144 0.00  0.00         2.98  0.007         6

Regression: Absolute scores

Model 0a: Disorder_hierarchical, Unadjusted

# A tibble: 8 x 6
  `Indepent variable`  Estimate `95% CI low` `95% CI up` SE    `p value`
  <fct>                <chr>    <chr>        <chr>       <chr> <chr>    
1 Intercept            2.25     " 1.42"      3.07        0.42  9.22e-08 
2 gad.sum_score_prepan 0.55     " 0.53"      0.56        0.01  0.00e+00 
3 Anxiety disorders    0.59     "-0.31"      1.49        0.46  0.20     
4 Eating disorders     0.79     "-0.98"      2.57        0.91  0.38     
5 OCDs                 1.42     " 0.19"      2.65        0.63  0.02     
6 PTSD                 1.15     "-0.27"      2.57        0.73  0.11     
7 Depressive disorders 1.51     " 0.68"      2.33        0.42  3.39e-04 
8 Bipolar disorders    1.86     " 0.98"      2.75        0.45  3.80e-05 

# A tibble: 1 x 6
  n     R2    `adjust. R2` F      `p value`    df
  <chr> <chr> <chr>        <chr>  <chr>     <dbl>
1 11144 0.32  0.32         743.84 0.00e+00      7

Model 0b: Overlapping disorders, Unadjusted

# A tibble: 14 x 6
   `Indepent variable`             Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                       " 2.98"  " 2.38"      3.58        0.30  1.34e-22 
 2 gad.sum_score_prepan            " 0.52"  " 0.50"      0.53        0.01  0.00e+00 
 3 Depressive and anxiety disorder " 0.83"  " 0.23"      1.42        0.30  0.006    
 4 Only depressive disorder        "-0.24"  "-0.88"      0.40        0.33  0.47     
 5 Only anxiety disorder           " 0.11"  "-0.56"      0.78        0.34  0.76     
 6 Eating disorders                " 0.73"  " 0.44"      1.02        0.15  6.17e-07 
 7 OCDs                            " 0.61"  " 0.34"      0.88        0.14  1.08e-05 
 8 Psychotic and bipolar disorder  "-0.53"  "-1.31"      0.26        0.40  0.19     
 9 Only psychotic disorder         " 0.39"  "-0.27"      1.05        0.34  0.24     
10 Only bipolar disorder           " 0.17"  "-0.23"      0.56        0.20  0.41     
11 PTSD                            " 0.75"  " 0.50"      1.01        0.13  5.96e-09 
12 ASD                             " 0.90"  " 0.38"      1.42        0.27  7.08e-04 
13 ADHD                            "-0.13"  "-0.84"      0.59        0.36  0.73     
14 Personality disorder            " 0.92"  " 0.56"      1.27        0.18  3.95e-07 

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 11165 0.01  0.01         11.04 2.61e-12      6

Model 1a: Disorder hierarchical, Adjusted

Model 0a + Gender_collapsed + age_category_collapsed_reordered + Ethnicity_collapsed + Time registration and baseline

# A tibble: 24 x 6
   `Indepent variable`     Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                   <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept               " 2.44"  " 1.48"      " 3.41"     0.49  7.57e-07 
 2 gad.sum_score_prepan    " 0.52"  " 0.50"      " 0.53"     0.01  0.00e+00 
 3 Anxiety disorders       " 0.56"  "-0.34"      " 1.46"     0.46  0.22     
 4 Eating disorders        " 0.67"  "-1.09"      " 2.44"     0.90  0.45     
 5 OCDs                    " 1.21"  "-0.01"      " 2.43"     0.62  0.05     
 6 PTSD                    " 1.42"  " 0.00"      " 2.84"     0.73  0.05     
 7 Depressive disorders    " 1.61"  " 0.79"      " 2.43"     0.42  1.29e-04 
 8 Bipolar disorders       " 2.08"  " 1.20"      " 2.97"     0.45  4.22e-06 
 9 Female                  " 0.52"  " 0.28"      " 0.76"     0.12  2.04e-05 
10 Non-binary/Self-defined " 0.81"  " 0.13"      " 1.50"     0.35  0.02     
# … with 14 more rows

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 11165 0.01  0.01         11.04 2.61e-12      6

Model 1b: Overlapping disorders, Adjusted

Model 0a + Gender_collapsed + age_category_collapsed_reordered + Ethnicity_collapsed + Time registration and baseline

# A tibble: 30 x 6
   `Indepent variable`             Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                       " 3.08"  " 2.30"      " 3.85"     0.40  1.02e-14 
 2 gad.sum_score_prepan            " 0.49"  " 0.48"      " 0.51"     0.01  0.00e+00 
 3 Depressive and anxiety disorder " 0.92"  " 0.33"      " 1.52"     0.30  0.002    
 4 Only depressive disorder        " 0.01"  "-0.63"      " 0.66"     0.33  0.96     
 5 Only anxiety disorder           " 0.07"  "-0.60"      " 0.74"     0.34  0.83     
 6 Eating disorders                " 0.61"  " 0.32"      " 0.89"     0.15  3.79e-05 
 7 OCDs                            " 0.53"  " 0.26"      " 0.80"     0.14  1.22e-04 
 8 Psychotic and bipolar disorder  "-0.26"  "-1.05"      " 0.52"     0.40  0.51     
 9 Only psychotic disorder         " 0.45"  "-0.20"      " 1.10"     0.33  0.18     
10 Only bipolar disorder           " 0.32"  "-0.07"      " 0.71"     0.20  0.11     
# … with 20 more rows

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 11165 0.01  0.01         11.04 2.61e-12      6

Model 2a: Disorder hierarchical, Adjusted with age²

Model 0a + Gender_collapsed + age² + Ethnicity_collapsed + Time registration and baseline

# A tibble: 18 x 6
   `Indepent variable`                   Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                                 <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                             "  2.18" "  1.24"     "  3.13"    0.48  6.32e-06 
 2 gad.sum_score_prepan                  "  0.52" "  0.50"     "  0.53"    0.01  0.00e+00 
 3 Anxiety disorders                     "  0.54" " -0.35"     "  1.44"    0.46  0.24     
 4 Eating disorders                      "  0.65" " -1.12"     "  2.41"    0.90  0.47     
 5 OCDs                                  "  1.22" "  0.00"     "  2.44"    0.62  0.05     
 6 PTSD                                  "  1.39" " -0.03"     "  2.81"    0.72  0.05     
 7 Depressive disorders                  "  1.58" "  0.76"     "  2.40"    0.42  1.69e-04 
 8 Bipolar disorders                     "  2.05" "  1.16"     "  2.94"    0.45  5.84e-06 
 9 Female                                "  0.50" "  0.26"     "  0.74"    0.12  3.76e-05 
10 Non-binary/Self-defined               "  0.79" "  0.11"     "  1.48"    0.35  0.02     
11 Age                                   "-71.59" "-81.58"     "-61.60"    5.10  2.00e-44 
12 Age^2                                 " -9.17" "-18.94"     "  0.60"    4.98  0.07     
13 Mixed or multiple ethnic origins      "  0.06" " -0.56"     "  0.69"    0.32  0.85     
14 Asian or Asian British                " -0.16" " -1.17"     "  0.86"    0.52  0.76     
15 African or African British            " -0.46" " -1.98"     "  1.06"    0.78  0.55     
16 Other ethnicity                       "  0.49" " -0.49"     "  1.46"    0.50  0.33     
17 Key worker                            " -0.34" " -0.53"     " -0.16"    0.10  3.12e-04 
18 Time difference: Lockdown to baseline "  0.00" " -0.06"     "  0.05"    0.03  0.95     

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 11165 0.01  0.01         11.04 2.61e-12      6

Model 2b: Overlapping disorders, Adjusted with age²

Model 0a + Gender_collapsed + age² + Ethnicity_collapsed + Time registration and baseline

# A tibble: 24 x 6
   `Indepent variable`             Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                       "  2.88" "  2.13"     "  3.63"    0.38  5.85e-14 
 2 gad.sum_score_prepan            "  0.49" "  0.48"     "  0.51"    0.01  0.00e+00 
 3 Depressive and anxiety disorder "  0.88" "  0.29"     "  1.47"    0.30  0.004    
 4 Only depressive disorder        "  0.00" " -0.64"     "  0.64"    0.33  0.99     
 5 Only anxiety disorder           "  0.04" " -0.63"     "  0.71"    0.34  0.90     
 6 Eating disorders                "  0.61" "  0.32"     "  0.90"    0.15  3.27e-05 
 7 OCDs                            "  0.52" "  0.25"     "  0.79"    0.14  1.44e-04 
 8 Psychotic and bipolar disorder  " -0.25" " -1.03"     "  0.54"    0.40  0.54     
 9 Only psychotic disorder         "  0.47" " -0.18"     "  1.12"    0.33  0.16     
10 Only bipolar disorder           "  0.32" " -0.07"     "  0.71"    0.20  0.11     
# … with 14 more rows

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 11165 0.01  0.01         11.04 2.61e-12      6

GAD prepan: sensitivity

Regression: Absolute scores

Model 1b: Overlapping disorders, Adjusted

Model 0a + Gender_collapsed + age_category_collapsed_reordered + Ethnicity_collapsed + Time registration and baseline

Subset of participants 3 months prior to pandemic

# A tibble: 30 x 6
   `Indepent variable`             Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                       " 2.34"  "-0.05"      4.72        1.21  0.05     
 2 gad.sum_score_prepan            " 0.54"  " 0.49"      0.60        0.03  1.22e-63 
 3 Depressive and anxiety disorder " 0.73"  "-0.92"      2.38        0.84  0.39     
 4 Only depressive disorder        "-0.27"  "-2.04"      1.51        0.90  0.77     
 5 Only anxiety disorder           " 0.11"  "-1.87"      2.08        1.01  0.92     
 6 Eating disorders                " 0.88"  "-0.16"      1.92        0.53  0.10     
 7 OCDs                            " 0.03"  "-1.02"      1.08        0.54  0.96     
 8 Psychotic and bipolar disorder  "-0.81"  "-3.33"      1.71        1.28  0.53     
 9 Only psychotic disorder         " 0.37"  "-2.23"      2.96        1.32  0.78     
10 Only bipolar disorder           "-0.53"  "-2.17"      1.11        0.84  0.53     
# … with 20 more rows

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 828   0.01  0.00         1.47  0.18          6

Subset of participants 6 months prior to pandemic

# A tibble: 30 x 6
   `Indepent variable`             Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                       " 2.12"  " 0.46"      3.79        0.85  0.01     
 2 gad.sum_score_prepan            " 0.54"  " 0.50"      0.58        0.02  1.01e-128
 3 Depressive and anxiety disorder " 0.82"  "-0.26"      1.90        0.55  0.14     
 4 Only depressive disorder        "-0.20"  "-1.39"      0.98        0.60  0.74     
 5 Only anxiety disorder           " 0.31"  "-0.99"      1.62        0.66  0.64     
 6 Eating disorders                " 0.55"  "-0.20"      1.30        0.38  0.15     
 7 OCDs                            " 0.58"  "-0.12"      1.28        0.36  0.11     
 8 Psychotic and bipolar disorder  "-0.47"  "-2.50"      1.56        1.03  0.65     
 9 Only psychotic disorder         " 0.47"  "-1.08"      2.02        0.79  0.55     
10 Only bipolar disorder           "-0.10"  "-1.15"      0.96        0.54  0.85     
# … with 20 more rows

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 1634  0.00  0.00         1.22  0.29          6

Subset of participants 9 months prior to pandemic

# A tibble: 30 x 6
   `Indepent variable`             Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                       " 2.69"  " 1.39"      " 4.00"     0.67  5.44e-05 
 2 gad.sum_score_prepan            " 0.53"  " 0.50"      " 0.57"     0.02  6.23e-209
 3 Depressive and anxiety disorder " 1.03"  " 0.18"      " 1.89"     0.44  0.02     
 4 Only depressive disorder        " 0.23"  "-0.70"      " 1.15"     0.47  0.63     
 5 Only anxiety disorder           " 0.50"  "-0.54"      " 1.54"     0.53  0.34     
 6 Eating disorders                " 0.91"  " 0.34"      " 1.49"     0.29  0.002    
 7 OCDs                            " 0.57"  " 0.03"      " 1.12"     0.28  0.04     
 8 Psychotic and bipolar disorder  "-0.68"  "-2.19"      " 0.82"     0.77  0.37     
 9 Only psychotic disorder         " 0.74"  "-0.47"      " 1.95"     0.62  0.23     
10 Only bipolar disorder           " 0.21"  "-0.59"      " 1.02"     0.41  0.60     
# … with 20 more rows

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 2710  0.01  0.00         2.99  0.007         6

##KT

GAD retro

Regression: Change scores

Model 0a: Disorder_hierarchical, Unadjusted

# A tibble: 7 x 6
  `Indepent variable`  Estimate `95% CI low` `95% CI up` SE    `p value`
  <fct>                <chr>    <chr>        <chr>       <chr> <chr>    
1 Intercept            2.58     " 2.44"      2.72        0.07  9.69e-272
2 Anxiety disorders    0.82     " 0.51"      1.13        0.16  1.80e-07 
3 Eating disorders     0.41     "-0.51"      1.33        0.47  0.38     
4 OCDs                 1.10     " 0.32"      1.88        0.40  0.006    
5 PTSD                 1.01     " 0.24"      1.78        0.39  0.01     
6 Depressive disorders 0.87     " 0.71"      1.04        0.09  4.00e-24 
7 Bipolar disorders    0.86     " 0.49"      1.24        0.19  5.71e-06 

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 18588 0.01  0.01         18.23 3.23e-21      6

Model 0b: Overlapping disorders, Unadjusted

# A tibble: 13 x 6
   `Indepent variable`             Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                       " 2.58"  " 2.45"      2.71        0.07  2.17e-303
 2 Depressive and anxiety disorder " 0.85"  " 0.68"      1.03        0.09  4.94e-21 
 3 Only depressive disorder        " 0.25"  " 0.01"      0.48        0.12  0.04     
 4 Only anxiety disorder           " 0.86"  " 0.58"      1.15        0.14  2.26e-09 
 5 Eating disorders                " 0.49"  " 0.22"      0.76        0.14  4.28e-04 
 6 OCDs                            " 0.05"  "-0.21"      0.31        0.13  0.72     
 7 Psychotic and bipolar disorder  "-0.20"  "-0.96"      0.56        0.39  0.61     
 8 Only psychotic disorder         "-0.28"  "-0.94"      0.38        0.34  0.41     
 9 Only bipolar disorder           "-0.21"  "-0.62"      0.19        0.21  0.30     
10 PTSD                            " 0.41"  " 0.17"      0.66        0.12  8.72e-04 
11 ASD                             " 0.36"  "-0.12"      0.83        0.24  0.14     
12 ADHD                            " 0.01"  "-0.63"      0.64        0.32  0.99     
13 Personality disorder            " 0.29"  "-0.06"      0.64        0.18  0.11     

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 18588 0.01  0.01         18.23 3.23e-21      6

Model 1a: Disorder hierarchical, Adjusted

Model 0a + Gender_collapsed + age_category_collapsed_reordered + Ethnicity_collapsed + Time registration and baseline

# A tibble: 22 x 6
   `Indepent variable`                   Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                                 <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                             " 4.16"  " 3.79"      " 4.52"     0.19  4.36e-106
 2 Anxiety disorders                     " 0.57"  " 0.25"      " 0.90"     0.17  5.51e-04 
 3 Eating disorders                      " 0.02"  "-0.91"      " 0.96"     0.48  0.96     
 4 OCDs                                  " 0.54"  "-0.25"      " 1.34"     0.41  0.18     
 5 PTSD                                  " 0.78"  "-0.02"      " 1.58"     0.41  0.06     
 6 Depressive disorders                  " 0.63"  " 0.45"      " 0.82"     0.10  3.17e-11 
 7 Bipolar disorders                     " 0.65"  " 0.26"      " 1.03"     0.20  9.77e-04 
 8 Female                                " 0.77"  " 0.57"      " 0.97"     0.10  2.26e-14 
 9 Non-binary/Self-defined               "-0.24"  "-0.90"      " 0.42"     0.34  0.48     
10 age_category_collapsed_reordered16-18 "-0.31"  "-0.74"      " 0.11"     0.22  0.15     
# … with 12 more rows

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 18588 0.01  0.01         18.23 3.23e-21      6

Model 1b: Overlapping disorders, Adjusted

Model 0a + Gender_collapsed + age_category_collapsed_reordered + Ethnicity_collapsed + Time registration and baseline

# A tibble: 29 x 6
   `Indepent variable`             Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                       " 3.94"  " 3.57"      " 4.31"     0.19  2.84e-95 
 2 Depressive and anxiety disorder " 0.65"  " 0.46"      " 0.85"     0.10  4.95e-11 
 3 Only depressive disorder        " 0.25"  " 0.00"      " 0.51"     0.13  0.05     
 4 Only anxiety disorder           " 0.63"  " 0.34"      " 0.93"     0.15  2.96e-05 
 5 Eating disorders                " 0.29"  " 0.01"      " 0.56"     0.14  0.04     
 6 OCDs                            "-0.09"  "-0.36"      " 0.17"     0.14  0.50     
 7 Psychotic and bipolar disorder  " 0.09"  "-0.67"      " 0.86"     0.39  0.81     
 8 Only psychotic disorder         "-0.05"  "-0.71"      " 0.62"     0.34  0.89     
 9 Only bipolar disorder           "-0.11"  "-0.52"      " 0.30"     0.21  0.60     
10 PTSD                            " 0.43"  " 0.18"      " 0.68"     0.13  6.93e-04 
# … with 19 more rows

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 18588 0.01  0.01         18.23 3.23e-21      6

Model 2a: Disorder hierarchical, Adjusted with age²

Model 0a + Gender_collapsed + age² + Ethnicity_collapsed + Time registration and baseline

# A tibble: 16 x 6
   `Indepent variable`                   Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                                 <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                             "  2.11" "  1.67"     "  2.55"    0.22  4.09e-21 
 2 Anxiety disorders                     "  0.54" "  0.13"     "  0.95"    0.21  0.010    
 3 Eating disorders                      "  0.38" " -0.80"     "  1.56"    0.60  0.52     
 4 OCDs                                  "  0.30" " -0.73"     "  1.34"    0.53  0.57     
 5 PTSD                                  "  1.17" "  0.06"     "  2.28"    0.57  0.04     
 6 Depressive disorders                  "  0.95" "  0.72"     "  1.18"    0.12  1.36e-15 
 7 Bipolar disorders                     "  1.00" "  0.56"     "  1.44"    0.22  7.90e-06 
 8 Female                                "  0.64" "  0.41"     "  0.86"    0.11  2.52e-08 
 9 Non-binary/Self-defined               " -0.60" " -1.34"     "  0.14"    0.38  0.11     
10 Age                                   "-78.52" "-90.65"     "-66.40"    6.19  1.10e-36 
11 Age^2                                 "  7.26" " -4.22"     " 18.74"    5.86  0.22     
12 Mixed or multiple ethnic origins      " -0.67" " -1.30"     " -0.04"    0.32  0.04     
13 Asian or Asian British                " -0.22" " -1.14"     "  0.69"    0.47  0.63     
14 African or African British            " -0.67" " -2.12"     "  0.79"    0.74  0.37     
15 Other ethnicity                       "  0.16" " -0.97"     "  1.29"    0.58  0.78     
16 Time difference: Lockdown to baseline " -0.07" " -0.12"     " -0.02"    0.03  0.007    

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 13420 0.01  0.01         30.14 4.04e-36      6

Model 2b: Overlapping disorders, Adjusted with age²

Model 0a + Gender_collapsed + age² + Ethnicity_collapsed + Time registration and baseline

# A tibble: 23 x 6
   `Indepent variable`             Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                       "  1.99" "  1.55"     "  2.43"    0.22  5.70e-19 
 2 Depressive and anxiety disorder "  0.95" "  0.71"     "  1.18"    0.12  5.37e-15 
 3 Only depressive disorder        "  0.48" "  0.17"     "  0.79"    0.16  0.002    
 4 Only anxiety disorder           "  0.60" "  0.22"     "  0.98"    0.19  0.002    
 5 Eating disorders                "  0.26" " -0.05"     "  0.58"    0.16  0.10     
 6 OCDs                            " -0.27" " -0.57"     "  0.03"    0.15  0.08     
 7 Psychotic and bipolar disorder  "  0.33" " -0.58"     "  1.23"    0.46  0.48     
 8 Only psychotic disorder         " -0.06" " -0.81"     "  0.68"    0.38  0.87     
 9 Only bipolar disorder           " -0.12" " -0.56"     "  0.33"    0.23  0.61     
10 PTSD                            "  0.58" "  0.30"     "  0.86"    0.14  5.60e-05 
# … with 13 more rows

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 13420 0.01  0.01         30.14 4.04e-36      6

Regression: Absolute scores

Model 0a: Disorder_hierarchical, Unadjusted

# A tibble: 8 x 6
  `Indepent variable`  Estimate `95% CI low` `95% CI up` SE    `p value`
  <fct>                <chr>    <chr>        <chr>       <chr> <chr>    
1 Intercept            3.53     3.40         3.67        0.07  0.00e+00 
2 gad.sum_score_retro  0.61     0.59         0.62        0.01  0.00e+00 
3 Anxiety disorders    1.81     1.52         2.11        0.15  1.59e-33 
4 Eating disorders     1.26     0.40         2.13        0.44  0.004    
5 OCDs                 2.70     1.97         3.44        0.38  5.91e-13 
6 PTSD                 2.16     1.43         2.88        0.37  5.60e-09 
7 Depressive disorders 2.28     2.11         2.45        0.09  1.24e-151
8 Bipolar disorders    2.94     2.57         3.30        0.18  7.53e-57 

# A tibble: 1 x 6
  n     R2    `adjust. R2` F       `p value`    df
  <chr> <chr> <chr>        <chr>   <chr>     <dbl>
1 18588 0.32  0.32         1276.84 0.00e+00      7

Model 0b: Overlapping disorders, Unadjusted

# A tibble: 14 x 6
   `Indepent variable`             Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                       3.75     " 3.62"      3.88        0.07  0.00e+00 
 2 gad.sum_score_retro             0.55     " 0.54"      0.57        0.01  0.00e+00 
 3 Depressive and anxiety disorder 2.31     " 2.14"      2.49        0.09  3.78e-148
 4 Only depressive disorder        0.67     " 0.45"      0.89        0.11  2.55e-09 
 5 Only anxiety disorder           1.89     " 1.63"      2.16        0.14  4.35e-44 
 6 Eating disorders                0.78     " 0.53"      1.03        0.13  1.83e-09 
 7 OCDs                            0.75     " 0.50"      0.99        0.12  2.15e-09 
 8 Psychotic and bipolar disorder  0.10     "-0.60"      0.81        0.36  0.77     
 9 Only psychotic disorder         0.03     "-0.58"      0.65        0.31  0.92     
10 Only bipolar disorder           0.16     "-0.21"      0.54        0.19  0.40     
11 PTSD                            0.84     " 0.61"      1.07        0.12  3.49e-13 
12 ASD                             1.02     " 0.58"      1.46        0.22  6.19e-06 
13 ADHD                            0.44     "-0.15"      1.03        0.30  0.14     
14 Personality disorder            1.09     " 0.76"      1.42        0.17  6.62e-11 

# A tibble: 1 x 6
  n     R2    `adjust. R2` F      `p value`    df
  <chr> <chr> <chr>        <chr>  <chr>     <dbl>
1 23100 0.15  0.15         660.52 0.00e+00      6

Model 1a: Disorder hierarchical, Adjusted

Model 0a + Gender_collapsed + age_category_collapsed_reordered + Ethnicity_collapsed + Time registration and baseline

# A tibble: 24 x 6
   `Indepent variable`     Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                   <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept               " 5.56"  " 5.21"      " 5.92"     0.18  3.96e-203
 2 gad.sum_score_retro     " 0.56"  " 0.54"      " 0.57"     0.01  0.00e+00 
 3 Anxiety disorders       " 1.48"  " 1.18"      " 1.79"     0.15  1.29e-21 
 4 Eating disorders        " 0.78"  "-0.09"      " 1.65"     0.44  0.08     
 5 OCDs                    " 1.91"  " 1.17"      " 2.65"     0.38  4.63e-07 
 6 PTSD                    " 1.98"  " 1.24"      " 2.73"     0.38  1.86e-07 
 7 Depressive disorders    " 2.07"  " 1.89"      " 2.25"     0.09  2.61e-109
 8 Bipolar disorders       " 2.77"  " 2.40"      " 3.14"     0.19  2.68e-49 
 9 Female                  " 0.76"  " 0.57"      " 0.94"     0.09  7.47e-16 
10 Non-binary/Self-defined " 0.28"  "-0.34"      " 0.90"     0.31  0.37     
# … with 14 more rows

# A tibble: 1 x 6
  n     R2    `adjust. R2` F      `p value`    df
  <chr> <chr> <chr>        <chr>  <chr>     <dbl>
1 23100 0.15  0.15         660.52 0.00e+00      6

Model 1b: Overlapping disorders, Adjusted

Model 0a + Gender_collapsed + age_category_collapsed_reordered + Ethnicity_collapsed + Time registration and baseline

# A tibble: 30 x 6
   `Indepent variable`             Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                       " 5.56"  " 5.22"      " 5.91"     0.18  3.62e-213
 2 gad.sum_score_retro             " 0.51"  " 0.50"      " 0.53"     0.01  0.00e+00 
 3 Depressive and anxiety disorder " 2.15"  " 1.96"      " 2.33"     0.09  2.61e-112
 4 Only depressive disorder        " 0.79"  " 0.55"      " 1.02"     0.12  3.46e-11 
 5 Only anxiety disorder           " 1.61"  " 1.33"      " 1.89"     0.14  2.57e-30 
 6 Eating disorders                " 0.52"  " 0.26"      " 0.77"     0.13  7.02e-05 
 7 OCDs                            " 0.57"  " 0.33"      " 0.82"     0.13  5.06e-06 
 8 Psychotic and bipolar disorder  " 0.42"  "-0.28"      " 1.12"     0.36  0.24     
 9 Only psychotic disorder         " 0.26"  "-0.36"      " 0.87"     0.31  0.41     
10 Only bipolar disorder           " 0.30"  "-0.08"      " 0.67"     0.19  0.13     
# … with 20 more rows

# A tibble: 1 x 6
  n     R2    `adjust. R2` F      `p value`    df
  <chr> <chr> <chr>        <chr>  <chr>     <dbl>
1 23100 0.15  0.15         660.52 0.00e+00      6

Model 2a: Disorder hierarchical, Adjusted with age²

Model 0a + Gender_collapsed + age² + Ethnicity_collapsed + Time registration and baseline

# A tibble: 18 x 6
   `Indepent variable`                   Estimate  `95% CI low` `95% CI up` SE    `p value`
   <fct>                                 <chr>     <chr>        <chr>       <chr> <chr>    
 1 Intercept                             "   2.86" "   2.45"    "   3.27"   0.21  3.52e-42 
 2 gad.sum_score_retro                   "   0.55" "   0.53"    "   0.57"   0.01  0.00e+00 
 3 Anxiety disorders                     "   1.53" "   1.14"    "   1.91"   0.19  4.59e-15 
 4 Eating disorders                      "   1.29" "   0.21"    "   2.37"   0.55  0.02     
 5 OCDs                                  "   2.09" "   1.14"    "   3.05"   0.49  1.75e-05 
 6 PTSD                                  "   2.58" "   1.56"    "   3.61"   0.52  7.56e-07 
 7 Depressive disorders                  "   2.53" "   2.31"    "   2.75"   0.11  1.56e-106
 8 Bipolar disorders                     "   3.23" "   2.81"    "   3.64"   0.21  5.28e-52 
 9 Female                                "   0.58" "   0.37"    "   0.79"   0.11  3.42e-08 
10 Non-binary/Self-defined               "  -0.04" "  -0.72"    "   0.64"   0.35  0.91     
11 Age                                   "-121.62" "-133.06"    "-110.18"   5.84  8.19e-95 
12 Age^2                                 "  10.80" "  -0.08"    "  21.68"   5.55  0.05     
13 Mixed or multiple ethnic origins      "  -0.45" "  -1.03"    "   0.13"   0.30  0.13     
14 Asian or Asian British                "   0.08" "  -0.76"    "   0.93"   0.43  0.85     
15 African or African British            "  -0.07" "  -1.41"    "   1.27"   0.68  0.92     
16 Other ethnicity                       "   0.55" "  -0.49"    "   1.59"   0.53  0.30     
17 Key worker                            "  -0.05" "  -0.22"    "   0.12"   0.09  0.57     
18 Time difference: Lockdown to baseline "  -0.02" "  -0.07"    "   0.02"   0.02  0.34     

# A tibble: 1 x 6
  n     R2    `adjust. R2` F      `p value`    df
  <chr> <chr> <chr>        <chr>  <chr>     <dbl>
1 16839 0.19  0.19         657.64 0.00e+00      6

Model 2b: Overlapping disorders, Adjusted with age²

Model 0a + Gender_collapsed + age² + Ethnicity_collapsed + Time registration and baseline

# A tibble: 24 x 6
   `Indepent variable`             Estimate  `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>     <chr>        <chr>       <chr> <chr>    
 1 Intercept                       "   2.95" "   2.55"    "  3.35"    0.20  7.90e-47 
 2 gad.sum_score_retro             "   0.51" "   0.50"    "  0.53"    0.01  0.00e+00 
 3 Depressive and anxiety disorder "   2.57" "   2.34"    "  2.79"    0.11  2.83e-109
 4 Only depressive disorder        "   1.14" "   0.85"    "  1.42"    0.14  4.22e-15 
 5 Only anxiety disorder           "   1.70" "   1.35"    "  2.05"    0.18  7.64e-22 
 6 Eating disorders                "   0.58" "   0.30"    "  0.87"    0.15  7.13e-05 
 7 OCDs                            "   0.40" "   0.12"    "  0.67"    0.14  0.004    
 8 Psychotic and bipolar disorder  "   0.39" "  -0.44"    "  1.21"    0.42  0.36     
 9 Only psychotic disorder         "   0.34" "  -0.34"    "  1.02"    0.35  0.32     
10 Only bipolar disorder           "   0.32" "  -0.09"    "  0.73"    0.21  0.12     
# … with 14 more rows

# A tibble: 1 x 6
  n     R2    `adjust. R2` F      `p value`    df
  <chr> <chr> <chr>        <chr>  <chr>     <dbl>
1 16839 0.19  0.19         657.64 0.00e+00      6

GAD retro: sensitivity

Regression: Change scores

Model 0a: Disorder_hierarchical, Unadjusted

# A tibble: 7 x 6
  `Indepent variable`  Estimate `95% CI low` `95% CI up` SE    `p value`
  <fct>                <chr>    <chr>        <chr>       <chr> <chr>    
1 Intercept            2.16     " 1.00"      3.32        0.59  2.65e-04 
2 Anxiety disorders    1.11     "-0.16"      2.38        0.65  0.09     
3 Eating disorders     1.99     "-0.40"      4.38        1.22  0.10     
4 OCDs                 1.16     "-0.57"      2.88        0.88  0.19     
5 PTSD                 1.80     "-0.09"      3.70        0.97  0.06     
6 Depressive disorders 1.19     " 0.03"      2.36        0.60  0.05     
7 Bipolar disorders    1.05     "-0.20"      2.30        0.64  0.10     

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 8058  0.00  0.00         0.94  0.47          6

Model 0b: Overlapping disorders, Unadjusted

# A tibble: 13 x 6
   `Indepent variable`             Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                       " 2.55"  " 1.72"      3.38        0.42  1.97e-09 
 2 Depressive and anxiety disorder " 0.73"  "-0.11"      1.57        0.43  0.09     
 3 Only depressive disorder        " 0.15"  "-0.76"      1.07        0.47  0.74     
 4 Only anxiety disorder           " 0.84"  "-0.11"      1.78        0.48  0.08     
 5 Eating disorders                " 0.55"  " 0.16"      0.93        0.20  0.006    
 6 OCDs                            "-0.21"  "-0.57"      0.15        0.18  0.25     
 7 Psychotic and bipolar disorder  "-0.30"  "-1.39"      0.79        0.56  0.59     
 8 Only psychotic disorder         "-0.16"  "-1.07"      0.74        0.46  0.73     
 9 Only bipolar disorder           "-0.29"  "-0.83"      0.24        0.27  0.28     
10 PTSD                            " 0.39"  " 0.06"      0.73        0.17  0.02     
11 ASD                             " 0.22"  "-0.49"      0.92        0.36  0.55     
12 ADHD                            " 0.12"  "-0.86"      1.09        0.50  0.81     
13 Personality disorder            " 0.42"  "-0.06"      0.90        0.24  0.09     

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 8058  0.00  0.00         0.94  0.47          6

Model 1a: Disorder hierarchical, Adjusted

Model 0a + Gender_collapsed + age_category_collapsed_reordered + Ethnicity_collapsed + Time registration and baseline

# A tibble: 22 x 6
   `Indepent variable`                   Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                                 <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                             " 2.52"  " 1.18"      " 3.86"     0.68  2.21e-04 
 2 Anxiety disorders                     " 0.99"  "-0.27"      " 2.25"     0.64  0.12     
 3 Eating disorders                      " 1.69"  "-0.69"      " 4.06"     1.21  0.16     
 4 OCDs                                  " 0.90"  "-0.81"      " 2.61"     0.87  0.30     
 5 PTSD                                  " 1.94"  " 0.06"      " 3.82"     0.96  0.04     
 6 Depressive disorders                  " 1.23"  " 0.06"      " 2.39"     0.59  0.04     
 7 Bipolar disorders                     " 1.17"  "-0.08"      " 2.41"     0.63  0.07     
 8 Female                                " 0.67"  " 0.33"      " 1.01"     0.17  1.22e-04 
 9 Non-binary/Self-defined               "-0.64"  "-1.56"      " 0.28"     0.47  0.17     
10 age_category_collapsed_reordered16-18 "-0.22"  "-0.99"      " 0.54"     0.39  0.57     
# … with 12 more rows

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 8058  0.00  0.00         0.94  0.47          6

Model 1b: Overlapping disorders, Adjusted

Model 0a + Gender_collapsed + age_category_collapsed_reordered + Ethnicity_collapsed + Time registration and baseline

# A tibble: 29 x 6
   `Indepent variable`             Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                       " 2.81"  " 1.73"      " 3.89"     0.55  3.49e-07 
 2 Depressive and anxiety disorder " 0.71"  "-0.14"      " 1.55"     0.43  0.10     
 3 Only depressive disorder        " 0.37"  "-0.54"      " 1.29"     0.47  0.43     
 4 Only anxiety disorder           " 0.67"  "-0.27"      " 1.62"     0.48  0.16     
 5 Eating disorders                " 0.37"  "-0.02"      " 0.76"     0.20  0.06     
 6 OCDs                            "-0.27"  "-0.63"      " 0.09"     0.18  0.15     
 7 Psychotic and bipolar disorder  " 0.03"  "-1.05"      " 1.12"     0.55  0.95     
 8 Only psychotic disorder         " 0.11"  "-0.80"      " 1.02"     0.46  0.81     
 9 Only bipolar disorder           "-0.17"  "-0.70"      " 0.37"     0.27  0.54     
10 PTSD                            " 0.46"  " 0.12"      " 0.80"     0.17  0.008    
# … with 19 more rows

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 8058  0.00  0.00         0.94  0.47          6

Regression: Absolute scores

Model 0a: Disorder_hierarchical, Unadjusted

# A tibble: 8 x 6
  `Indepent variable`  Estimate `95% CI low` `95% CI up` SE    `p value`
  <fct>                <chr>    <chr>        <chr>       <chr> <chr>    
1 Intercept            4.62     3.56         5.69        0.54  2.04e-17 
2 gad.sum_score_retro  0.52     0.50         0.55        0.01  0.00e+00 
3 Anxiety disorders    1.15     0.00         2.31        0.59  0.05     
4 Eating disorders     2.74     0.56         4.93        1.11  0.01     
5 OCDs                 2.21     0.64         3.78        0.80  0.006    
6 PTSD                 2.11     0.39         3.84        0.88  0.02     
7 Depressive disorders 1.73     0.66         2.79        0.54  0.001    
8 Bipolar disorders    2.17     1.03         3.31        0.58  1.86e-04 

# A tibble: 1 x 6
  n     R2    `adjust. R2` F      `p value`    df
  <chr> <chr> <chr>        <chr>  <chr>     <dbl>
1 8058  0.20  0.20         290.49 0.00e+00      7

Model 0b: Overlapping disorders, Unadjusted

# A tibble: 14 x 6
   `Indepent variable`             Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                           <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept                       " 5.45"  " 4.69"      6.21        0.39  2.82e-44 
 2 gad.sum_score_retro             " 0.48"  " 0.46"      0.51        0.01  0.00e+00 
 3 Depressive and anxiety disorder " 0.89"  " 0.13"      1.65        0.39  0.02     
 4 Only depressive disorder        "-0.64"  "-1.46"      0.18        0.42  0.13     
 5 Only anxiety disorder           " 0.67"  "-0.18"      1.52        0.43  0.12     
 6 Eating disorders                " 0.88"  " 0.53"      1.23        0.18  8.46e-07 
 7 OCDs                            " 0.51"  " 0.19"      0.84        0.17  0.002    
 8 Psychotic and bipolar disorder  "-0.33"  "-1.31"      0.65        0.50  0.51     
 9 Only psychotic disorder         " 0.14"  "-0.68"      0.95        0.42  0.74     
10 Only bipolar disorder           " 0.07"  "-0.41"      0.55        0.25  0.77     
11 PTSD                            " 0.84"  " 0.54"      1.15        0.16  6.51e-08 
12 ASD                             " 0.86"  " 0.23"      1.49        0.32  0.008    
13 ADHD                            " 0.52"  "-0.36"      1.39        0.45  0.25     
14 Personality disorder            " 1.26"  " 0.83"      1.69        0.22  1.31e-08 

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 9688  0.01  0.01         11.18 1.82e-12      6

Model 1a: Disorder hierarchical, Adjusted

Model 0a + Gender_collapsed + age_category_collapsed_reordered + Ethnicity_collapsed + Time registration and baseline

# A tibble: 24 x 6
   `Indepent variable`     Estimate `95% CI low` `95% CI up` SE    `p value`
   <fct>                   <chr>    <chr>        <chr>       <chr> <chr>    
 1 Intercept               " 4.78"  " 3.55"      " 6.01"     0.63  2.77e-14 
 2 gad.sum_score_retro     " 0.50"  " 0.48"      " 0.53"     0.01  0.00e+00 
 3 Anxiety disorders       " 1.16"  " 0.02"      " 2.31"     0.58  0.05     
 4 Eating disorders        " 2.59"  " 0.43"      " 4.74"     1.10  0.02     
 5 OCDs                    " 2.00"  " 0.44"      " 3.55"     0.79  0.01     
 6 PTSD                    " 2.51"  " 0.80"      " 4.22"     0.87  0.004    
 7 Depressive disorders    " 1.99"  " 0.93"      " 3.05"     0.54  2.34e-04 
 8 Bipolar disorders       " 2.54"  " 1.40"      " 3.67"     0.58  1.17e-05 
 9 Female                  " 0.53"  " 0.22"      " 0.84"     0.16  7.40e-04 
10 Non-binary/Self-defined " 0.01"  "-0.82"      " 0.84"     0.42  0.99     
# … with 14 more rows

# A tibble: 1 x 6
  n     R2    `adjust. R2` F     `p value`    df
  <chr> <chr> <chr>        <chr> <chr>     <dbl>
1 9688  0.01  0.01         11.18 1.82e-12      6

Model 1b: Overlapping disorders, Adjusted

Model 0a + Gender_collapsed + age_category_collapsed_reordered + Ethnicity_collapsed + Time registration and baseline